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Flow past a flat plate at low Reynolds numbers 


By E. JANSSEN 


General Electric Company, San Jose, California 


(Received 26 July 1957) 


SUMMARY 

‘The flow past a flat plate at Reynolds numbers in the range 
0-1 to 10-0 is investigated by an analogue method. ‘The solution 
gives the stream function and the vorticity in the flow field 
surrounding the plate. From these are obtained the local coefficient 
of friction, the pressure distribution along the plate, and the total 
drag coefficient. ‘The drag coefficient approaches the analytical 
values of Haaser (1950) and of ‘Tomotika & Aoi (1953) as the 
Reynolds number decreases toward 0-1. ‘The drag coefficient 
approaches the Blasius solution as the Reynolds number increases. 
At Reynolds number 10-0 the drag coefficient is still above the 
Blasius value, but is below the value obtained experimentally by 
Janour (1951). The difference from the experimental result is 
attributed for the most part to truncation error. 


INTRODUCTION 

The Blasius solution of the Prandtl boundary layer equation, for the 
case of flow past a thin flat plate, is valid only for a boundary layer whose 
thickness is very small relative to distance from the leading edge. ‘The 
results of the Blasius analysis are therefore invalid at low Reynolds numbers. 
Kuo (1953) has investigated the flow past a flat plate for Reynolds numbers 
as low as 10 by a perturbation method. ‘Tomotika & Aoi (1953), Imai 
(1954), and others have investigated the flow past a flat plate for Reynolds 
numbers up to about 1-0, by linearizing the flow equation in accordance 
with the Oseen approximation. Haaser (1950) has computed the skin 
friction distribution for this configuration, using Carrier’s (1953) modifica- 
tion of the Oseen linearization for Reynolds numbers up to 5-0. The 
object of this paper is to investigate, by an analogue method, the flow past 
a finite flat plate-at Reynolds numbers in the range 0-1 to 10-0, using the 
exact non-linear formulation of the problem, and to compare the results 
with those of the foregoing and with experiment. 

The equation governing the two-dimensional motion of an incom- 
pressible fluid with constant viscosity is 


Vy Vo = yV2C 


where ¢ is the scalar vorticity, given by 
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The defining equations of a stream function # are 
Ors Ors 
“= - ay? v= =? (3) 
In Cartesian coordinates, these equations may be written 


ren ee 2 


ox dy dy Ox 


V4 = C. (5) 
Equations (3) satisfy the continuity equation which, therefore, does not 
need to be expressed explicitly. 

Equations (4) and (5) have the form of the Poisson equation. ‘These 
particular equations have been solved for the case of flow past a circular 
cylinder with zero circulation at Reynolds numbers in the range 1 to 1000, 
by numerical techniques (Thom 1933; Kawaguti 1953; and Allen & 
Southwell 1955). The procedure outlined in this paper gives a method of 
solution by means of two resistance networks which is applicable to an 
extensive class of problems in viscous flow. It will be applied here 
specifically to the case of flow past a flat plate. 


Transformation of equations 

It is convenient at the outset to render the equation dimensionless, and 
to consider transformation from the physical plane to some other plane. 
Let the physical plane be the z-plane, and let the velocity of the undisturbed 
flow be parallel to the x-axis. ‘Then 


, 


2’ = x +t’ =1 


ft pt0 


e 
Here the prime notation signifies dimensional quantities. ‘The hydro- 
dynamic equations, as well as the expression for z, are rendered dimensionless 
by making the following substitutions : 











e u’ v p’ 
2=5 u=—- o= = ——, 
ry uw ca 
py’ af i ae N at. %o 
Ve oe a yf = aera c= 3 2 = ’ 
4 “7%, u", v 


where 7, is some characteristic length, e.g. a cylinder radius, u, is the 
undisturbed velocity parallel to the x-axis, v’ is the kinematic viscosity 
and N» is the Reynolds number. 

The transformation t = f(z) can be employed to map the field on to the 
t-plane where 


t = €+1». (6) 
If the transformation is conformal, equations (4) and (5) become 

es Nr | Cus ara ous 0b) si 

Vil = ae 5 Se ae (7) 


and 
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A2 a2 
oO 


“ 
= ArT 5 

0& Oo 
and q, the transformation modulus, is the ratio of a length in the ¢-plane 
to the corresponding length in the z-plane, and is a function of position, 
that is 


where > 
Vi 





dt 
dz 


q = = g(x,y) = q(g,%). 








Solution of the Poisson equation by a network analogue 
The general form of the Poisson equation 
V*d = f(x,y) (9) 
may be expressed approximately in the difference form 
= 1 72 , y 
bo = bu —4@ Ff (X03 o)s (10) 
where the field has been divided into square meshes of length a on a side 
figure 1). is the value of the function in question at a particular mesh 
0 1 
point, and ¢,, is the average for the nearest four neighbouring mesh points. 

















a 
$2 
| o3 $o $) | 
| 
a 
4 { 
—a— Cl —e 





Figure 1. Representation of field by a net. 


Now consider a resistance network having a geometrically similar mesh 
pattern (figure 2), with resistances R; between adjacent mesh points (nodes) 
as shown. In addition to the current flowing through these resistors, 
current flows into each of the nodes from an external source. ‘The sum of 
the currents flowing into a particular node must be zero from Kirchhoff’s 
law of node currents: 

Co— 41 - &— fn , hg —h3 Og % 
R, R, R R, 
dR; = R;—R, 
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where dR; is the deviation of the jth resistor with respect to the first resistor, 
The analogy that exists between the function ¢ of the system described by 
equation (10) and the voltage e of the network is apparent. Ifascale factor K 
is chosen so that e = K¢ at corresponding points, then 7 can be adjusted so 
that 
4 2 i § "qe 
page Soy ee haath (12) 
ee j 1jaz ANG 1 
Note that when the resistor deviation is zero the first two terms on the right 
of equation (12), which are in fact error terms due to this deviation, drop 
out. By the method of this paper, however, it is not necessary that the 
resistor deviation be zero, as the current i, itself is not measured, and the 
network resistor deviation contributes nothing to the computational error. 





Figure 2. Resistance network pattern geometrically similar to net of figure 1. 


The Poisson equations whose solutions are desired are equations (7) 
and (8). The difference forms of these equations are, respectively, 


bo = Sur — ae Nn {(ts — bs)(S2— Sa) — (He — a(S — Ss)} (13) 
and e. ; 
in by = by < c. (14) 
where 
Cu = (4) + Co + C3 + C4) by = (py ae iy Tr ibs + 4). 


There is a network for each equation. ‘The voltages at the boundaries 
of the two networks are fixed to correspond to the boundary conditions 
of the problem. Current must then be supplied to each anode in the 
s-network proportional (except for the resistor deviation terms of (12)) 
to ¢ for that point, and to each node in the ¢-network proportional to 
(ab, — th3)( Co — C4) — (Wo — y)(C, — %g)} for that point. ‘To accomplish this 
quantities, called the ‘residuals’ in relaxation procedures, are defined by 


9 


q » 
KR, : td 4h 2 Lo (15) 








} 
| 
| 
} 
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for the J-net, and by 
R. at 40 ii 4¢o = sNets = bs)(Co = C4) ai (i. i. Ba)(Or 7 Cs)5 (16) 


for the ¢-net. ‘The current to each node in the %-network is adjusted 
manually until #, goes to zero, and the current to each node in the ¢-network 
is adjusted manually until #, goes to zero. Z#, and Z, may be computed 
electronically by using components which change sign, add or multiply 
(the theory is clearly presented in a number of texts, for example see 
Korn & Korn (1952) or Soroka (1954)). A schematic diagram of the networks 
and associated voltage sources, current controls, etc., is shown in figure 3. 
One arrangement of electronic components for computing residuals is 
shown schematically in figure 4. 
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Figure 3. Schematic diagram of analogue. 


The solution to equations (13) and (14) is obtained after a relaxation 
procedure involving the adjustment of the currents to each of the nodes in 
the two networks in several successive steps. ‘This is the approximate 
solution to equations (7) and (8), the degree of approximation improving 
with the fineness of mesh. A mesh size should be used such that any 
decrease in the size will produce no change in the solution within the 
desired limits of accuracy. 


Boundary conditions (general) 


In the case of a flat plate located at or near the origin in an otherwise 
uniform flow field, % is essentially at its uniform field value, and ¢ is zero, 
far away from the origin. ‘Thus, this part of the boundary, henceforth 
referred to as the outer boundary, is chosen far away and this may be taken 
care of by making the network quite large. It may also be handled by using 
a transformation which permits a reasonably small network to be used as 
was done in obtaining the results reported in this paper. 
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ys will be known at every place on the boundary and ¢ will be known at 
every place except where the boundary is solid (that is, along the flat plate), 
Consider a solid boundary which is parallel to the €-axis (figure 5). The 
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Figure 4. Assembly of electronic components for computing residuals. 
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Figure 5. Solid boundary parallel to €-axis. 


value of ¢ at the point 0 on the boundary is desired. Points 1 and 3 are 
on the boundary and point 2 is a short distance inside the boundary at a 
position on a normal through point 0. ‘The expression for ¢ in terms of 
follows from equation (8): 


: de of (ozs O7ubs\ 7 e 
{= @lVi¥), = ail (2), + (=) |. (17) 


A difference form for equation (17) may be obtained by using a Taylor 
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series expansion about point 0. ‘The expression for the stream function 


at point 2 is 

oy a* (d7ys a® (By a (org F 
Yo = Wot a TS wee ae Tae oe T Aa aa + O(a). 
Yo (; Se x5; ), 6 (5), 24\ 0n*/o ore 


From equation (8), 
Oy 1 
==) +—&, 
og" 0 qi : 


~ - au aclg 
(),- Let Le Ca)e 
oy (OM) (AAG) _ (Ul 
on} 0 - 04 0 On? 1? -  0g2 0 


Along the solid boundary % is constant and u is zero. Hence, all derivatives 
of y with respect to €, and the first derivative of % with respect to 7, are zero. 
i. is now given by 


ce Co : a® (0(C/q") a o*(C/q°) o°(¢/q") 
Ye = si Po + 2 @ at on 0 i 24 “—e'), 7 Oe 0 ; — 


0 























Neglecting all terms past the second on the right (that is, all terms of order a® 
or higher), and rearranging, one obtains 
Co = Af (5 — ty). (19) 
Thom (1933) and Allen & Southwell (1955) have used equation (19), 
and their results check well with experiment. It was possible for them to 
decrease the mesh constant a in regions where (c¢/en), was large, thus 
keeping the third (and following) term on the right of equation (18) small 
and presumably negligible. Woods (1954), on the other hand, takes the 
third term into account in the following manner. ‘The Taylor series is 
again employed, this time to obtain an expression for the vorticity at 





point Z: 
ix, “46 O(f/q (C/¢ 
-= = -af = 4) + a (ote) 4 ’) +0 (a?). 
qs 45 on Jo 2\ on? /0 


Rearrangement and multiplication by ja? gives 


= q)\ @f @% a ‘¢ *(¢/9") 5 
6\ on Jy 6 q5 6 qi ie on? Jo nen 


and substitution for the third term in (18) then gives 


ac a* ¢ a c 
Loa = +> 4 ~ a eS | 5 20 
2 Por 3 ge 6 q5 mal t Al 0(a’). (20) 





Neglect of all terms in (20) of order a‘ or higher, and rearrangement, 
finally gives 


ae. SY ee tbe (21a) 


a“ 
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or, more generally, 

qe 
242, 
where ¢, is the vorticity at any point s on the solid boundary, ¢, is the stream 
function at that point, ¢,, is the iaiewns at a point p inside the boundary 
at a distance m ona narnia! through s, and y%, is the stream function at 
point p. 

Equation (21) is preferable to (19) for use with networks unless the mesh 
constant has been reduced in regions of the flow field where the vorticity 
gradient is steep. Equation (21) was used in obtaining the results reported 
in this paper. 


Si 


hus (21 b) 


s 


o] » 
Iq; 
—= (bs — ad. 
m2 (%,, Pe 


Values of ¢, for every point on the solid boundary are first assumed, 
then new values are determined after each relaxation of the two networks 
If succeeding values of ¢, are based on equation (21b), that 1s, 





u~ 


grey = SEE (yon _y, 2¢ C0), 
where the superscript in parentheses refers to the number of the step in 
yal iteration, the process will diverge. ‘To effect convergence a plot of 
{€.+(g2/2q5)6, }% vs bY? may be prepared and the graph of equation (21 b) 
superpose d. The intersection of the two curves will give the final value of 
{f.+(q?/2q7)¢,}. A supplementary plot of ¢) vs (¢.+(g5/2g;)G, } may 
assist in estimating the correct final value of ¢ 


Flow past a flat plate 
Consider the case of a finite flat plate (of infinitesimal thickness) in an 
infinite flow field, and parallel to the otherwise undisturbed velocity, with 
leading edge at (—1,0) and trailing edge at (+1,0). Let the trans- 
formation be 
z= coat, (22) 
for which the modulus is 
g = [sinh? €+ sin? y]-??. (23) 
This transformation has the following features: 
(1) the modulus is very small far away from the origin, thus permitting 
the transformation plane to be relatively small; 
(2) the modulus is very large in the region of the leading edge, thus 
literally magnifying the region of interest. 
Because the flow will be symmetrical only the half plane need be used. 
The transformation is illustrated in figure 6. 
The boundary conditions are: 


é = 0, 4s = 0, == 0; 
SS, y = —sinhé, siny, C= 0. 


y = U0, 7, ob = (), f= 0. 
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The solution obtained from the analogue represents % and ¢ as functions 
of position in the flow field (including the solid boundary). ‘The velocity 
is given by (3). ‘The equation for the pressure is 


b Ps 7 Os 2 6é | ous 2 of 
t , C ‘ , OY 
ee oe ( ee FN ae 2 
sue | y 3 € N R°) 2 } s ie 1] N R 0€ ~ (24) 


The shear stress is given by 























Zz 
itera 25) 
=> Rt (25) 
Z- PLANE +-PLANE 
id 
FLAT 
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FLAT 
PLATE 
_ 10) 
| foo . r@) -* 


Figure 6. ‘Transformation for a flat plate. 


The local coefficient of friction c, is equal numerically to twice the dimension- 
less shear stress of equation (25). Let the defining equation for the drag 
coefficient C, be 

D' = Cphl'p'ut, 
where D’ is the drag on both sides of the plate per unit span, and /’ is the 
dimension of the plate parallel to the flow. ‘Then Cp is given by 
1 0) 


| c;dx = —5| csiny dy. (26) 
1 / 2 


Cp = 


bol = 
hol 


Results and discussion 

Two networks were constructed, with outer boundary at €,, = ?7* and 
with mesh constant equal to {vf (a transition to $a was effected beyond 
&=17). The REAC of the Rand Corporation, Santa Monica, California, 
was made available to the author for the purpose of computing residuals, 
its components being connected as indicated in figure 4. Data obtained 
from the analogue in the form of potentials were readily converted to stream 
function and vorticity by the application of the appropriate scale factors. 

* Location of the outer boundary was based on a consideration of the distance 
from the plate at which the maximum vorticity in a laminar wake would have become 


small enough to neglect. 
+ Considerations based on limiting the truncation error in ¢ led to the value 


ia for the mesh constant. ‘This was increased by a factor of 5 for reasons of economy. 
Any error present in the solution is presumed to be largely truncation error. 
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Solutions were obtained for Reynolds numbers 0-1, 1-0, and 10-0. The 
corresponding values of ys and ¢ are given in tables 1 and 2. ‘The differenta- 
tion and integration involved in determining the pressure and drag coefficient 
were done numerically. 

The local coefficient of friction is plotted vs x, for the three Reynolds 
numbers, in figure 7. c, as given by the Blasius solution is also plotted. ‘The 
analogue solution approaches the Blasius solution as the Reynolds number 
increases. ‘The influence of the trailing edge, significant at very low 
Reynolds numbers, may be seen distinctly in the analogue solution. ‘The 
Blasius solution, of course, does not take the effect of the trailing edge 
into account. 

















































































































































































































I | ED 9 i as & i | } i i 
ANALOG | ANALOG T] 
BLASUS =e — BLASUS — ——— — 
>—+4 }_4 
75 4 24 
50 \ / 
, Np 20. VA | 
25 16 
0 oe oe Ce ee en \ NR "0.1 
Pl | 
"10 
R 
8 
I 
me v4 NR #100 
\ ses Ia 
P-Pw o | \ a sae eee eee See ee 
8 \ N iu an — — =— 
Nelo — ~\ 
4r v 
Sch 
me ee bee am lee anlae oe == an 
0 — = =o 
— 
7 N -16 
aly aL 
N. =10.0 
me x ame re -24 
‘\ ~ 
ed ee + mee ee ee 
-1.0 - 0.5 re) 0.5 10 -10 -0.5 fe) 0.5 1.0 
x x 
Figure 7. Local friction coefficient vs x. Figure 8. Pressure vs x. 


‘The pressure at the surface of the plate was found by first integrating the 
terms of (24) from the outer boundary to the plate along the path » = 37 
(and thereby determining the pressure at the mid-point of the plate), and 
then integrating along the surface of the plate (that is, along the path € = 0). 
The pressure is plotted vs x, for the three Reynolds numbers, in figure 8. 
The value given by the Blasius solution is zero everywhere. It may again 
be noted that the analogue solution approaches the Blasius solution as the 
Reynolds number increases. 
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The drag coefficient Cp is given by (26). One should know c¢; at every 
point along the x-axis from x = —1 to x= +1 in order to evaluate the 
integral. But the analogue gives values only from x = —0-924 to 
x = +0-924. The fact that c; is not known in the regions adjacent to the 
leading and trailing edges is a serious deficiency. However, some corre- 
lation is established between the solutions at the three Reynolds numbers 
by plotting c, vs the ‘local’ Reynolds numbers Nz, and Nyy, in 














n Nr | = 0 dor hor Sar dar ar 7 oa 
O-t | 0 0-022 0-08 0-21 — — — — 

ta 1-0 0 0-030 0-11 0:28 — — — - 
10-0 0 0-076 0-27 0-57 — — — — 
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Sar — 0 0-050 0-20 0-50 1-04 3°21 8-04 17-94 
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oar — 0 0-060 0-25 0-61 = = a 
= 0 0-131 0-54 1°21 = co 














—~ 0 0-049 0-2 25°3 
dir — | 0 0-060 0-25 0-62 1-30 4-14 11-10 25-3 
— 0 0-116 0-4 253 








= 0 0-045 0-18 0-46 -- = = 
7 — 0 0-051 0-20 0-51 = = 
— 0 0-079 0-33 0-81 — — — = 








— 0 0-035 0-14 0:34 0-71 2°31 6:63 17-94 
ls — 0 0-036 0-14 0-35 0-70 | —2:20| —6:-27 17-94 
— 0 0-042 0-17 0-41 0-79 3-64 8°38 17:94 











— 0 0-020 0-07 0-19 — — = = 
Lar — 0 0-021 0-07 0-18 — = == = 
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Table 1. Stream function. 


figures 9 and 10 respectively. Np, is the local Reynolds number based 
on distance from the leading edge and N,,_,) is the local number based 
on distance from the trailing edge. On the basis of these plots a straight 
line extrapolation (on logarithmic coordinates) is employed to obtain 
values of c, in the regions —1 < « < —0-924 and +0-924 <x < +1, 
the slope of the straight line being taken as — 0-577. 
Cp is now given by 

1 0°924 1 cl 

Cp =5 | c;dx—5| c,sinn dy+ 5 | cy dx. 


l ga ~ 0-924 


+47 


bole 
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The first term on the right is the contribution of the 3-8°%, segment of the 
surface adjacent to the leading edge, the last term is the contribution of the 
3-8, segment adjacent to the trailing edge, and the middle term is the 
contribution of the 92:4°%, of the surface which lies between. 
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Let these 
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Table 2. Vorticity. 
Nr c DI c ‘D2 ( ‘D3 C 'D 
| 0-1 7-09 30-6 6:76 44-45 
| 1-0 1:07 3-78 0-75 5-60 
| 10-0 0:28 0:79 0:07 1-14 
| 
| | 
Table 3. Drag coefficient. 


contributions to C, be denoted respectively by Cp,, Cp3, and Cops, so that 
Cp = Cor + Coe + Cos: 
is given in table 3. 
According to this analysis the contribution of the leading edge segment 
increases from 16%, to 25°, of the total drag, the contribution of the trailing 


The drag coefficient, together with its three parts, 
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edge segment decreases from 15°, to 6%, and the contribution of the 
surface between stays nearly constant at 68°, or 69°%% as the Reynolds 
number is increased from 0-1 to 10-0. This is in contrast to the Blasius 






















































































100 
a i) es 
~ CANIN tees 
jones 
| 4 Np* 0:1 BLASIUS 
~~ — 
~~ aati 
| ~~ ™, 
Su, Mie 
O | ~ 
~ a Paid Np =1.0 
‘| — i" 
q 
a ” ~d 
~ _— 
> Np 710.0 
~ 
> — 
~ 
~ 
| 
ooo 00! 0.1 1.0 10.0 
N 
Re 
Figure 9. Local coefficient of friction vs N,,. 
00 
10 
‘ 
0.00! ool 01 
Nace~ XO 


Figure 10. Local coefficient of friction vs Np; 


solution according to which the contribution of the same leading edge 
segment is constant at 13-75% of the total drag and the contribution of the 
trailing edge segment is constant at 1-35°% of the total drag, for all Reynolds 


numbers. 
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The drag data of table 3 are compared with the experimental values 
reported by Janour (1951) in figure 11. Drag data reported by Sherman 
(1952) are also shown although these were obtained under compressible 
flow conditions at Mach numbers of approximately 0-2 and 0-6. Also shown 
are the drag coefficient for the Blasius solution and the drag coefficients 
as determined by Tomotika & Aoi (1953) and Haaser (1950), the latter two 
being based respectively on the Oseen approximation and on Carrier’s 
(1953) modification of the Oseen approximation. ‘The drag coefficient as 
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Figure 11. Drag coefficient vs Reynolds number. 


determined from the analogue solution approaches the Blasius solution 
as the Reynolds number increases, and approaches the solutions of Haaser 
and of ‘Tomotika & Aoi as the Reynolds number decreases. In the range 
where comparison with experiment is possible, the drag coefficient lies 
below the values of Janour and Sherman. ‘This difference from experiment 
is believed to be due chiefly to truncation error. 


‘he author wishes to thank the Rand Corporation, Santa Monica, 
California, for making available to him their analogue computing facilities. 
The solutions reported here were all obtained at Rand. ‘The resistance 
networks and associated switching equipment were made possible by the 
University of California Research Grant No. 1115. 
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turbulent boundary layer 
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SUMMARY 

This paper describes the results of further experimental 
investigation of the turbulent boundary layer with zero pressure 
gradient. Measurements of autocorrelation and of space-time 
double correlation have been made respectively with single 
hot-wires and with two hot-wires with the separation vector in 
any direction. Space-time correlations reach a maximum for 
some optimum delay. In the case of two points set on a line 
orthogonal to the plate, the optimum delay 7; is not zero. In 
the general case it is equal to the corresponding delay 7,, increased 
by compensating delay for translation with the mean flow. ‘Taylor’s 
hypothesis may be applied to the boundary layer at distances 
from the wall greater than 3°, of the layer thickness. Space-time 
isocorrelation surfaces obtained with optimum delay have a large 
aspect ratio in the mean flow direction, even if they are relative to 
a point close to the wall (0-035); the correlations along the mean 
flow then retain high values on account of the large scale of the 
turbulence. 


1. INTRODUCTION 

These experiments, following our previous investigations (Favre, 
Gaviglio & Dumas 1952, 1953a, b, 1954c, 1957), deal with longitudinal 
components of turbulent velocity in the boundary layer on a flat plate, 
with zero pressure gradient, set in a low turbulence wind tunnel. ‘The 
plate is made of a glass 210 cm in length, 80 cm in width, with a profiled 
leading edge and prolonged with an adjustable flap. ‘Transition is obtained 
by roughness elements beginning at the leading edge. 

The velocity outside the boundary layer is 12-20 msec, and the 
conventional boundary-layer thickness 6 is 17-5,33 and 34 mm for distances 
from the leading edge of 79, 185 and 194 cm (Reynolds numbers 
R, = 14500, 27900 and 28 000). 

Measurements have been made with hot-wire anemometers, and the 


apparatus needed for time correlations, autocorrelations, and space-time 
double correlation coefficients (see Favre 1946, 1948; Favre, Gaviglio & 
Dumas 1953a). ‘lhe space-time correlation coefficient R,,(7,X,, Xo, X3) 
for a point A is measured by setting one hot-wire at fixed position A and 
another at any other position B; the two signals are recorded simultaneously 
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and played back with a relative delay 7). ‘The autocorrelation coefficient 
R(T, 0, 0,0) is obtained when the points A and B are the same (a single 
hot-wire). ’ and y will be used to denote the distances from the wall of 
wires A and B respectively. X, is the component of their separation 
parallel to the plate (positive in the direction of the mean velocity); 
X,; = y’—y; Xz, is the separation orthogonal to X, and X3. 

The influence on the measurements of the disturbances produced by 
the upstream probe was checked, and lessened by improvement of probes. 
With regard to space-time correlations the disturbance reached 13° when 
the two hot-wires were set on the same mean streamline and close 
by one another (X,/5 = 0-25 for y’/5 = 0-77). No disturbance appeared 
if distance X,/5 was greater than 0-4 for y'/5 = 0-03, and 1-5 for 

= 0:77; this is not a certain criterion for disturbances to be negligible, 
as it may be accepted as such, at least if the upstream wire is a little more 
distant from the wall than the downstream wire. 

‘The mean velocity profiles are in agreement with the accepted shapes, 
and the universal profile proposed by Clauser (1954). 

2. SPACE-TIME CORRELATIONS 

The space-time correlation coefficient R,,(7,X,, Xz, X3) with respect to 
a point A fixed by its distance y’ to the wall and its distance 2’ to the 
leading edge, is said to be measured: 

‘longitudinally’ when X, = X; = 0 (X, and T variable) ; 
‘transversely’ when X, = X, = 0 (X, and T variable) ; 
‘longitudinally and transversely’ when X,=0 (Xj, X3 and 7 
variable) ; 
‘longitudinally, transversely and laterally’ when X,, X., X3 and 7 
are variable. 
Numerous measurements were made; figure 1 represents as an example 
some values of space-time correlations measured ie eae transversely 
and laterally with regard to a fixed point A (y’ = 1mm; 2: 185 cm). It 
is found that the correlation R,, changes with dalees T and reaches a maximum 
for an optimum value of delay to be denoted by 7’, in general and by 7’; in 


the particular case of transverse correlations. 


Optimum delay T, for transverse space-time correlation R,,(T, 0,0, X3) 

Figures 2 and 3 represent (V,,—V,,)T;/6 (V,, and V, are mean velocities 
at the positions A and Bb) sana when the trannition is obtained by 
preturbulence (Favre, Gaviglio & Dumas 1954a) and by roughness 
respectively; the differences are small. Comparison at distances s = 79cm 
and z = 194 cm (the difference being about 45 times the mean thickness 
of boundary layer) shows that for the same values of y’/6 and y/6 the delay 7; 
is proportional to 6. It will be noticed that the variation of 7; as a function 
of z is very small by comparison with the time required for movement 
with the mean flow over the same range of z. Figure 4 shows values of 7; 
for other distances y’/d of fixed wire from the wall. 
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Figure 2. Optimum delay T; for transverse space-time correlations in the boundary; 
layer on a flat plate (transition by preturbulence). 
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Figure 3. Optimum delay 7; for transverse space-time correlations in the boundary; 
layer on a flat plate (transition by roughness). 
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Optimum delay T,,, for transverse and longitudinal space-time correlations 
Ry(T,X,, 0, Xs) 

As previously (Favre, Gaviglio & Dumas 1954b, 1957), we verified 
in these investigations (especially at distances X/5 = 4) that, as a first 
approximation, the space-time correlation of velocity fluctuations at the 
fixed point A and at another point B reaches a maximum for a delay T,, 
equal to the ‘initial’ optimum delay 7; (which gives maximum transverse 
space-time correlation for a point C set on the mean streamline passing 











Figure 4. Optimum delay 7; for transverse space-time correlation in the boundary 


layer on a flat plate (transition by roughness). Z 185 cm. 
curve : a b ci 
vy /d 0-03 0-15 0:77 


through the point B, with X, = 0) increased by a delay 7), to compensate 
for the translation due to the mean flow; i.e. 
sue 
Fo See Foes where 7, = es : 

ViatVep 
Chen, when we pick out several positions B’, B” for point B on the same 
mean streamline, values of 7°, differ from each other by an amount 
. ( CB’ Cr’ + ae 
~\Va + Veq V4 2 “Vat Ve 
since streamlines are close to lines of equal velocity. 

We have published (Favre, Gavigllio & Dumas 
measured and computed from the preceding relation for y’/d = 0-06 in 
the case of transition obtained by preturbulence, and for y’/é = 0-24 in 


, 


1957) values of 7, 


the case of transition obtained by roughness elements. In figure 5 are given 
values of 7, corresponding to y’/d5 equal to 0-03, 0-15 and 0-77. ‘The above 
relation is found to be satisfactory to a first approximation when the upstream 
wire is a little more distant from the wail than the downstream wire* 
(X,X, > 0). Yet this relation is lacking in accuracy when the upstream 
wire is closer to the wall than the downstream wire (X,X, < 0). This 


may be due to the disturbances produced by the upstream probe. 
; : if I 


* Except in the case corresponding to y’/d = 0:03, X, = 76:2 mm, for which 


the accuracy of measurement of 7’, is poor. 
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(c) 


Optimum delays T,,, for longitudinal and transverse space-time correlations. 


0-77. 


- 0-15, (c) y’/d 
© T,, computed. 


0-03, (b) y’/5 


(a) y'/d 


0, X, variable. 
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Autocorrelation R,,(T, 0, 0, 0) and longitudinal space correlations 
Ry, (0, X4, 0, 0); Taylor's (1938) hypothesis 

The previous experiments (Favre, Gaviglio & Dumas 1957) have 
shown ‘Taylor’s hypothesis to be in very satisfactory agreement with 
autocorrelation curves and longitudinal space correlation curves for the 
turbulence downstream of a grid of mesh size WM when the transformation 
T = X,/V is used, at least up to distances beyond which the correlation 
becomes too small for a conclusive comparison (X,/M > 1:5). Ananalogous 
comparison has been made for points inside the boundary layer at several 
distances from the wall, viz. y’/5 = 0-77, 0-24*, 0-15 and 0-03 (see figure 6). 
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Figure 6. Autocorrelations and longitudinal space correlations; ''aylor’s hypothesis. 
autocorrelations, ® longitudinal space correlations. 


‘The longitudinal space correlation curves (the streamlines are almost 
parallel to the plate) are seen to lie close to the autocorrelation curves 
for y’/d >0-03 when the transformation T = X,/V,,isused. ‘Uhe differences 
noticed, negligible in the central part of boundary layer (v'/5 = 0-24), have 
opposite signs on the two sides of this zone. 

It seems that to a first approximation ‘laylor’s hypothesis may be applied 
to the turbulent boundary layer, for distances from the wall greater than 


about 6°, of the boundary layer thickness. 


Space-time correlations with optimum delay 
We have given in a previous paper (Favre, Gaviglio & Dumas 1957) 
longitudinal and transverse space-time isocorrelation curves R,,(T,,,,X,,0,X3) 


7 


In this case transition is obtained by preturbulence, but the intermittency 
factor is negligible. 
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measured with optimum delay T.,, at distances y’/6 of the fixed wire A from 


ss the wall equal to 0-06, 0-24 and 0-71, and also at y’/6 > 1 in the case of 
measurements outside the boundary layer downstream of a grid producing 

- the preturbulence. 

h When X;, is given, the space-time correlation R,,(T7.,,,X4,0,X 3) with 

a optimum delay 7), reaches a maximum maximorum for one value of Xz. 


mm 


. Corresponding values of X, and X, determine the ‘maximum space-time 


, line with optimum delay’ or ‘maximum correlation line’ that goes through 
: the fixed point A. The present experiments were made in the case of 
l 


transition obtained by roughness and extended to three space dimensions, 
longitudinally, transversely, and laterally. 

In figures 7, 8, 9 we show space-time isocorrelation surfaces with 
optimum delay for positions of the fixed point A defined by y’/d = 0-03, 
0-15 and 0-77 respectively (z = 185 cm). ‘The symmetry with respect 
to X, having been checked, we may define as above the line of maximum 
correlation corresponding to each position of the fixed point A. Figure 10 
gives the space-time correlation coefficients with optimum delay 77, along 
the lines of maximum correlations, for positions of the fixed point A defined 
by y’/d = 0-03, 0-15 and 0-77. ‘These measurements were made with 
band-passes extending from 1 to 2500 Hz, and from 1 to 275 Hz by cutting 
off high frequencies. 

Space-time isocorrelation surfaces with optimum delay relating to a 
fixed point 4 have a great aspect ratio along the mean flow upstream and 
downstream, even at a small distance of the fixed point A from the wall 
(y’/6 = 0-03). ‘They are symmetrical with respect to X,, and their dimensions 
are of the same order transversely as laterally. ‘The lines of maximum 
correlation outside the boundary layer are the same as the mean streamlines, 
but they differ markedly inside the boundary layer. 





| 





| | | 
y'/d | 4 | 0-77 | 0-71 | 0-24 | 0-15 | 0-06 | 0-03 

| | | 

= rl 

‘Transition by preturbulence | 7:5 | | | | | 

R, 14.000 | to | —8-5| 2-4 | fm 
(band-pass 1—2500 Hz) 12 | 
| 


lransition by roughness 
2 7 | | ; 
R, 27 900 ei 2% 0-65 

H | 

(band-pass, 1—2500 Hz) | | 


‘Transition by roughness 
R, 27 900 10 3-6 1-6 
(band-pass, 1-275 Hz) 





Table 1. Values of X,/6 at which R,, = 0°5. 
‘The correlation coefficient with optimum delay retains high values for 
distances Y,/5 along the mean flow upstream and downstream which are 


appreciable at y’/5 > 0-03 and increase with y’/6. As an example, the 
value R,, = 0-5 is preserved as far as the distances indicated in table 1. 
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‘The distances are increased markedly by cutting off frequencies higher 
than 275 Hz, a result similar to that concerning the turbulence behind a 
grid (Favre, Gaviglio & Dumas 1954c). 
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Figure 10. Space-time correlations in the boundary layer on a flat plate. 
© Along the line of maximum correlation, with optimum delay 7’,,, band-pass 
1—2500 Hz. 
Along the line of maximum correlation, with optimum delay 7,,,, band-pass 
1-275 Hz. 


© Along a mean streamline, with zero delay. 


3. CONCLUSIONS 

The space-time correlation reaches a maximum for an optimum delay 7° 
When several alternative positions B’, B’ on a mean streamline are chosen for 
one of the points (B), the values of 7,, differ from each other by a 
compensating delay due to the motion from B’ to B” with the mean flow, 
the effective velocity of translation being taken as the average of the mean 
flow velocity at the point B and at the other point A. 

The values of the ‘initial’ optimum delay (corresponding to the case 


m* 


of the two points being set orthogonally to the plate) noted at two stations 
for which the distances from the leading edge are very different (the 
distances of the points from the wall being the same fraction of 6) are 
proportional to the boundary layer thickness 6. Variation of 7. with z 
is very small compared with the times corresponding to translation with 
the mean flow over the same range of x. ‘This seems to be due to the 
renewal of turbulence between the two distant positions. 

Taylor’s hypothesis may be applied, as a first approximation, to the 
boundary layer, for distances from the wall greater than about 6°, of the 
thickness 6. 

Space-time isocorrelation surfaces, with optimum delay relative to a 
fixed point A, have a great aspect ratio along the mean flow, upstream and 
downstream, even for small distances of the fixed point from the wall 
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(3% of 5). ‘Their dimensions are of the same order transversely and 
laterally. ‘lhe maximum correlation lines and mean streamlines are 
different inside the boundary layer, but are coincident outside. 

The correlation coefficient with optimum delay retains high values 
for distances upstream and downstream which are significant for y’/6 = 3% 
and increase in proportion to the distance of this fixed point from the wall. 
These distances are increased by cutting off the high frequency components. 

Along the mean streamlines the intensity of turbulence is practically 
constant and the correlation with optimum delay shows, on the one hand, 
inasmuch as it retains high values, the long persistence of large scale 
turbulence translated by the mean motion, and, on the other hand, in the 
light of its diminution, the evolution of the turbulence and its renewal 
inside the boundary layer. 


‘These researches were made at the Laboratoire di) Mécanique de 
l’Atmosphere for the Office National d’Etudes et de Recherches Aéro- 
nautiques (O.N.E.R.A.) with the aid of the Ministére de |’Air, and of the 
Centre National de la Recherche Scientifique. ‘The authors are glad to 
acknowledge Professor L. S. G. Kovasznay’s advice in 1956. 
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1. INTRODUCTION AND SUMMARY 


In the following it will be shown that a simple argument based on the 
} 


use of the energy integral equation of the laminar boundary layer permits 
the derivation of a heat transfer formula valid for non-uniform temperature 
distribution and non-zero pressure gradients. ‘lhe formula is then shown 
to be identical in structure with Lighthill’s (1950) well-known results. 
Lighthill obtained his formula by solving the boundary layer equations in 
the von Mises form using operational methods. An elegant way to obtain 
the same results using exact similarity consideration was given by Lagerstrom 
(not yet published). ‘lhe derivation given here is probably the most simple- 
minded one and the method may be useful for other applications as well. 
Furthermore, it is shown that the approach can be slightly modified to 
permit application of the formula to flow near separation. ‘The latter 
result is applied to the Falkner—Skan solution for just separating flows and 
is found to ‘be in excellent agreement with the exact solutions. 
2. LAMINAR FLOW FAR FROM SEPARATION 
‘The energy integral equation relating the heat flux g,, at the surface 
to the enthalpy distribution /(y) through the boundary layer can be 
written 
(a) 
2 | pu(h—h.,) dy = —q,, (1) 
dx J 0 


where 6, p, u are boundary layer thickness, density, and x-component of 
the velocity vector, respectively. Heat produced by dissipation has been 
neglected. 

If the Prandtl number Pr is not very small*, g,, is determined by the 
velocity distribution in the immediate neighbourhood of the wall. ‘This is 
a well-known result first used by Fage & Falkner (1931). ‘Thus 

u = (0u/Oy) eY = (Trl bw) (2) 
where 7,, and yw, are shearing stress and viscosity at the wall. We will 
now express y in formula (2) in terms of / and q in the neighbourhood of 
the wall and then employ a change of variables in (1) such that g is considered 
as a function of #. ‘This transformation is patterned after Crocco’s method 
of using 7(x,u) instead of z(x,y). We write 

h=h,,+(ch/ey) wy = hy —(Cp/R)w Ie Y> 


*'The opposite case Pr <1 is even easier since one can replace pu by ps U in (1). 
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and, on combining this with (2), we have near the wall 


t= Pr, (h,, —h)/ay- (3) 


The heat flux g is given by 


=>] => (4) 


The integration in (1) is performed at constant x and we can replace dy 
by dh using (4). ‘Thus (1) can be written 
1 d (7, (** (h,,—h)\(h—h,)pk 
G 


ey og ; 
Pr dx |qu Jn, qc, 





dh ‘ = dw 
P| 
or, if we introduce the similarity variable 
h,,—h 


= 


ete 
where h,, stands for the enthalpy at the edge of the layer, we have 
| @ {totut Oe a ye \yx(1—n) , ) a 
ar Bs 5 TwhwPw (h,,. -h.)* | a fi dy f= hoe (5) 
Perdx| ¢ Pwlw/ 9/Fw 


We now make the similarity assumptions that py/(p,,j,,) and g/q,, are 
functions of 7 only, that is 


fkhi—-h 
q(x; h) = qulx)f\ 7 - } 


While this cannot be exact in all cases, it is probably very closely satisfied 





for a large range of conditions, a fact which will be borne out by a comparison 
of the final results with known solutions. If the similarity is granted then 
the integral in (5) is a constant x, say, to be evaluated later. The result (5) 
becomes a very simple differential equation from which q,, as a function 


of 7,. and of Ah = h,,—h,, is easily obtained. We have from (5) 
ld (Ty bw Pw AA*\ Pr 
q si d \ \ ted } % . 
Multiplication with 4/(7,,.,,,p,, 4h?) makes the left side a perfect differential, 


and hence we find that 


— — } oe i | Te ie ge | a p,,, Ah®) dx : 
q \=%/ L Jo a 
.@) 
2a \' ), 2/3 13 3 A L3 y | Lis \ 
] |} PY ahr ip ARON CT, Mi Py ENA) | (6) 
3 | | 
jf ont 


The result (6) permits the computation of q,, if 7,,(x) and 7,,(x) or h,,(x) 
are given 
‘To determine the constant « we have to choose a reasonable distribution 


Of GPyMi/Fw pe =f(n)*. However, since only the cube root of « appears, 


One can also use the known solution for a flat plate with uniform temperature 
to obtain a, 
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the results are very insensitive to the choice of the function f(y). If we 
base our approach again on Crocco’s work for the shearing stress we choose 


f(a) = V— 7"), 
a simple function with the proper behaviour at the limits. ‘lhe constant « 
then has the value 1-47 = 0-215. ‘To compare (6) with Lighthill’s (1950) 
formula it is easiest first to find q,,(x) from (6) if Ah vanishes in the range 
(0,€) and then jumps to a constant value Ah, for greater values of x. In 
this way (6) can be written in the form of a Stieltjes integral like Lighthill’s 
formula, and the comparison is as follows. 
Lighthill’s result (1950, equation (29), with pu assumed constant) is 


g-13 pPr 1/3 “a t 3 J: 
rat Gu(%) = "ra (oe) \ tT (Xx) } | | \ iT w()$ ae | dT, (€). 


3 Ne J 
Equation (6) transforms into exactly the same form but with a numerical 
constant (3x)? instead of 9-13/2!. With the value « = 0-215 determined 
above we obtain 0-524 for the constant compared with Lighthill’s value 
of 0-539. ‘The agreement is thus better than 3°... Even with the worst 
possible choice of f(y), namely f(y) = 1, our constant would be reduced to 
only 0-490. 
3. LAMINAR FLOW NEAR SEPARATION 

It is now interesting to modify the approach to handle the flow near 
the separation point where 7,.-> 0. Here (2) obviously fails to represent u 
and the approximation becomes bad. However, it is an easy matter to try 
an approximation which represents uw near a separation point. Here, 


u == Ly?(d2u/dy?), (7) 
but p,(07u/dy*),,, = dp/dx. (8) 
Hence u (3 . 2p ) dp dx, 


and, proceeding as before, we obtain 


/ dp J T- dp 13 “| —1/4 
\14 Al a,% 4 | < J 44 | \ 
— 7, = (ep) Pi ( u p,, —- Ah Ue Dens Ah4 dx | F (Y) 
: ax ax | 
where f stands for the integral 
tT pe nh EF 
—— ——_——- dn, 
» Pu qq 
[he result (9) can be compared with the Falkner—Skan solution for just- 
separating flow. Lighthill gives a table showing the ratio (Nusselt 
number)/4 (Reynolds number) for Falkner Skan flow and Pr = 0-7. The 
value given for this ratio for separating Falkner—Skan flow is 


Nu//(Re,,) = 0-438. 


lhe value obtained from (9), again using f(y) = +/(1— 7”) in the evaluation 


ot 6, comes out to be 


Nu/4/(Re,.) = 0-448, 
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0 


and agrees also within better than 3°. In addition (9) shows explicitly 
the dependence of heat transfer on the Prandtl number near separation, 


4, EXTENSION TO HIGH MACH NUMBER FLOW 
The exact form of equation (1), valid for arbitrary Mach number, is 
d (° ES ‘ani 
Th | [h+3u?—(h, + $U*)]pu dy Gy (10) 
dx . r : 


Introducing the recovery enthalpy 4,, we may write the integral in the 


U 


form 
| [(h—h,)+(h,—h,,) — }(U? —u*)]pu dy. 


It is often possible to neglect the terms 
h h l( U2 u”) 


and extend the derivations to the high Mach number case simply by replacing 
h,. by h,. It is also possible of course to use the full equation (10), but 
then the differential equation becomes more difficult and the essential 


simplicity of the approach is lost. 
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Turbulent flow in a stably stratified atmosphere 
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SUMMARY 

Fluctuations of velocity and temperature which occur in 
a turbulent flow in a stably-stratified atmosphere far from 
restraining boundaries are discussed using the equations for 
the turbulent intensity and for the mean square temperature 
Huctuation. From these, an equation is derived for the flux 
Richardson number in terms of the ordinary Richardson number 
and some non-dimensional ratios connected with the turbulent 
motion. It 1s shown that the interaction between the temperature 
and velocity fields imposes on the flux Richardson number an 
upper limit of 0-5, and on the ordinary Richardson number a 
limit of about 0-08. If these values are exceeded, no equilibrium 
value of the turbulent intensity can exist and a collapse of the 
turbulent motion would occur. Although the analysis applies 
strictly only to a homogeneous non-developing flow, it should 
have approximate validity for effectively homogeneous, developing 
Hows, and the predictions are compared with some recent 


observations of these flows. 


1. INTRODUCTION 


‘The study of turbulent flow in a stably-stratitied fluid may be helped 
by the recognition of some different types of irregular flow which may 
occur. If the direct influence of the buoyancy forces on the motion is 
small, the motion will resemble ‘ordinary’ turbulence and will be 
characterized by high rates of energy dissipation, and of momentum and 
mass transport. If the buoyancy forces are dominant, the motion may be 
an irregular collection of gravity waves with low rates of energy dissipation 
and transport. In fows unrestrained by boundaries, both types of motion 
are found, characteristically turbulent motion near the origin of the flow 
and wave-like motion far downstream. ‘There is here a gradual transition 
from the region of turbulent flow to the region of wave flow, which makes 
dificult the definition of the limits of either flow: nevertheless, the two 
Hows are so distinct that no common description is likely to be valid. 

In this paper, the characteristically turbulent flow in a stably-stratified 
fluid far from restraining boundaries is discussed, using ideas and 
generalizations taken over from our knowledge of turbulent flow of constant 
density, with the purpose of obtaining a criterion for the continued 
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existence of fully turbulent flow. Naturally, the results can only be valid 
while the flow is fully turbulent, and the accuracy of the analysis depends 
on the occurrence of a breakdown in the mechanism of ‘ordinary’ 
turbulent flow before buoyancy forces become dominant. 

The methods used are by no means new, and might be described as 
mutton (mixing-length theory), dressed as lamb by recasting it in terms 
of similarity hypotheses. ‘The approach is somewhat similar to that 
used by Ellison (1957) who considered a boundary layer transmitting 
constant shear stress and constant heat flux. ‘This is an example of the 
class of turbulent flows which have complete homogeneity in the direction 
of flow allied to essential inhomogeneity in the direction of shear. All 
such flows are restrained by the fluid boundaries, while the unrestrained 
flows considered in this paper are substantially homogeneous in the 
direction of shear and inhomogeneous in the direction of flow. 

Although this paper is concerned primarily with flows in which radiative 
exchange of heat is negligible, the possibility of radiative transfer has been 
kept in mind and its effects will be considered in more detail in a following 
paper. 


2. CONDITIONS FOR SIMILARITY OF THE TURBULENT MOTION 


The attempt to find the criterion which determines whether turbulent 
flow does or does not occur in a particular flow system is a part of the 
general problem of establishing the conditions under which flows with 
geometrically similar boundary conditions are dynamically similar, the 
independent non-dimensional parameters whose equality permits dynamical 
similarity losing one degree of freedom when they specify a marginally 
stable flow. With the earth’s atmosphere in mind, we consider only flows of 
a nearly perfect gas, moving with velocity variations small compared with 
the local speed of sound, with temperature variations small compared withthe 
absolute temperature, and with a length scale small compared with 
the scale height of the atmosphere. If effects of the earth’s rotation 
are negligible, it may be shown (Batchelor 1953) that dynamical similarity 
of flows with geometrically similar boundary conditions, defined quanti- 
tatively by the three scales up, /), 4, referring to velocity variation, length, 
and potential temperature respectively, depends on equality of the three 
parameters tg /)/v (Reynolds number), g69 /y/( 7) u5) (Richardson number), 
v/« (Prandtl number), and equality or irrelevance of parameters depending 
on the radiative characteristics of the gas. Here vis the kinematic viscosity, 
« the thermometric conductivity, and 7) the mean absolute temperature. 
For flows of essentially the same gas the Prandtl number does not vary, 
and experience of shear flows of high Reynolds numbers shows that 
viscosity has no effect on the large-scale components of the motion which 
contain nearly all the energy. It follows that similarity depends only on 
the scale Richardson number and on the radiation parameters. 


No further simplification can be made by dimensional analysis alone 
without considering the dynamics of the flow. ‘lo the approximation 
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implied by the restrictions set out above, the Navier-Stokes equations of 


motion are D(U, +u;) : Sip 1 
Dt T p 

and the continuity equation is 
a(u,+ U,) 

Ox, 
where U; and u; are the 7-components of the mean fluid velocity and its 
fluctuation, @ is the temperature fluctuation about the mean 7, P and p 
are the mean, and the fluctuation about the mean, of the pressure difference 
from the ‘hydrostatic’ pressure P, defined by 0P,)/dx,; = —pg;, —g; is 
the acceleration vector of the gravitational field, p is the mean density. 
To this approximation, the direct effects of density variations on the 
velocity field are completely represented by the bouyancy term g,6/T, 
although the distribution of the buoyancy forces will depend on the 
interaction of the turbulent motion and the heat sources of the flow. 
The action of these buoyancy forces may be described in several ways but 
a fundamental effect is the addition or subtraction of energy from the 
turbulent motion. ‘The equation for the turbulent kinetic energy per unit 
mass 1s 
aU 4 W192 

& (ig) + um, = + me (1% + «pu. ) + U, ee = a bu; + vu;V?u;, (2.3) 


€ 


.) 


o( P +; 
OP? + (Ui +), (2.1) 


= 0, (2.2) 








where g? = u;u; and repeated suffices indicate summation over the three 
possible values of the repeated suffix. On the left-hand side of this equation, 
the first term is zero if the flow is statistically steady and the third and fourth 
terms represent transport of mechanical energy from one part of the flow 
to another, repectively transport by diffusive movements and through 
advection by the mean flow. Diffusive transport is usually effective only 
in the direction at right angles to the mean flow and can be eliminated 
from the equation by integrating over a whole section of the flow. Energy 
transport by advection cannot be eliminated in this way and it is of 
importance in all developing flows, e.g. boundary layers, jets and mixing 
layers, but its effects cannot usefully be considered without specifying 
the flow and it will be neglected in the following discussion. For this 
reason the results obtained apply strictly only to non-developing flows. 

Three terms now remain in equation (2.3), they represent energy 
production by transfer from the mean flow and energy loss through work 
done against buoyancy forces and through viscous dissipation. ‘To the 
extent that energy sources control its nature, the flow depends on the 
relative magnitudes of these three terms, or the ratio of the total loss of 
energy by buoyancy forces, 

aa 
| pou dA, 
to the total production by transfer from the mean flow, 
| wu; U, pn dA, 


Ox 


~ é 
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where the integrals extend over a whole section of the flow. This ratio, 


[ Ry], es 





2 } “—0U, ; 
pou; dA /| U;U, 7m 4A, (2.4) 


is the flux Richardson number for the whole section of the flow. ‘The 
ordinary flux form of the Richardson number, 





a ou. 
R, = 7 bu, || « u, = b (2.5) 


is the relevant /ocal parameter only if turbulent transport of energy can be 
neglected. ‘This analysis has been set out here to show that [R;], describes 
the relative effects of inertial and buoyancy forces on the turbulent motion, 
whether or not radiative transfer is an important element in the problem. 

Since radiative transfer plays no direct part in the dynamics of the 
turbulent motion or in the turbulent energy balance, it is a plausible 
assumption that the flux Richardson number will describe the motion with 
or without radiative transfer of heat. It does not follow that the condition 
for the maintenance of turbulent motion must be expressible as a single 
critical value of the flux Richardson number. Obviously the flux Richardson 
number cannot exceed one (turbulent dissipation of energy cannot change 
sign) but laboratory observations of flows unaffected by radiative transfer 
indicate a critical value of the flux Richardson number substantially less 
than one. ‘The critical condition appears to arise from a failure to achieve 
equilibrium in the interactions between the temperature field and the 
turbulent motion, which may occur in the following way. In a flow with 
given gradients of mean temperature and velocity, production of turbulent 
energy is proportional to the Reynolds stress and approximately proportional 
to the turbulent intensity. Loss of energy through buoyancy forces is 
proportional to the vertical heat flux or, in the absence of radiation, to the 
square-root of the turbulent intensity. ‘Turbulent dissipation of energy 
is proportional to the }-power of the intensity. Clearly the ratio of the 
total energy loss to the energy production, considered as a function of the 
turbulent intensity, has a minimum value, and if this is greater than one no 
equilibrium intensity exists and the turbulence must decay. Introduction 
of radiative heat transfer into the system alters the dependence of 
turbulent heat transport on turbulent intensity and will lead to a change in 
the critical condition. Although the motion itself may be still determined 
by the flux Richardson number, the range of possible Richardson numbers 
will depend on the intensity of the radiative transfer. ‘The next sections 
are concerned with the quantitative expression of these notions. 
3. "THE TURBULENT HEAT TRANSFER 

‘To the approximation being used, the equation for the temperature of 
a fluid element is 

D(T'+ 9) +(U,+u)# = kVAT+0)+—, (3.1) 


Dt oc. Ply 
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where & is the net rate of heat gain by radiative processes. If this equation 
is multiplied by the temperature fluctuation and the mean value taken, an 
equation for the intensity of the temperature fluctuations is obtained, 








Pe ee . a) O48) 8 i 
—(10")+u (o : =) U, —— + ——(3@u,) = «0V?0— Be. (3.2) 
x OX, 
In this equation, the radiation term -#@/(pc,,) has been written as — B6, 
where B may be regarded as the logarithmic rate of cooling of a fluid element 
by radiation alone. This radiation term is in general a function of the whole 
temperature field and its exact value will be considered in a second paper. 
From equation (3.2), an estimate of the heat transport by convective 
movements can be obtained. From this point on, it is convenient to consider 
a stationary flow with vertical gradients of the horizontal mean velocity and 
of temperature. In addition, the transport terms will be omitted from the 
equations for the turbulent energy and for the intensity of the temperature 
fluctuations, which means that advection is neglected and that the quantities 
in the equations are suitably defined mean values tor the whole flow. Using 
the ordinary notation with the Oz axis vertical and mean velocities along 


the Ow axis, equation (3.2) for 54° becomes 
wey if | go = (w 12 = 
wi (5 4 ge) _ pee, (3.3) 
Os Cy EB. 
where the conduction term «0V?@ has been replaced by — $67(w*)!*/L,. 


This substitution is to some extent formal, but it expresses also the 
experimentally established result that the large-scale properties of turbulent 
flows are independent of the magnitudes of the viscosity and conductivity 
of the fluid. In constant density flows, the length L, is nearly equal to the 
integral scale of the turbulent motion. From this equation, 


OT /oz +g/c, 





P2\U2 — p (gy2\ni2 | 2 A) 
(0?) k,(w") TRENT, 7? (3.4) 
5(w*) / sot f 
where 
k, = \w8| /[w?6?]!?, 
and the heat transfer by turbulent convection is 
Aw =— |UW Fon a = (3.5) 
R,, \(w?)l2/L,+4+ 8 


where Rk, = uw w*. The flux Richardson number, defined by equation (2.5), 


1S 
ke 
Res 





(3.6) 


An expression for the flux Richardson number may be obtained in a 
form involving only gradients of mean values and some non-dimensional 
ratios. ‘This is done by using the equation for the kinetic energy of the 
velocity fluctuations which is, omitting energy transport terms in the same 
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way as transport terms were omitted in equation (3.3), 


~~ g— (w?)32 
—uw = = — —0w+ —— ae 
a 7 3S (3-7) 
where L, is the dissipation length scale and is nearly equal to the integral 
scale in constant density flows. Combining this equation with (3.5) leads 
to a quadratic equation for the turbulent heat transport, whose solution 
may be written in the non-dimensional form 


(g/T)bw Ow | ( Lok; ib) ‘| 
R = 4H| 1-(1-12— a 3.8 
1 ua BU jaz Lk? H? (3.8) 
where +i, §£ R 
Halt Sy soc (3.9) 


and i g oT | g au . 
rior  & oz 


The number H is a measure of the ratio of the logarithmic rate of radiative 
cooling of a fluid element to the mean rate of shear. In the absence of 
radiative effects, it equals one. 


4. CONDITIONS FOR THE MAINTENANCE OF TURBULENT MOTION 

In the previous section, an equation was obtained relating the turbulent 
heat transport to the gradients of mean velocity and potential temperature, 
to the radiative properties of the fluid, and to some non-dimensional ratios 
describing the turbulent motion. ‘The validity of conclusions drawn from 
this equation depends only on further assumptions about the variation of 
the non-dimensional ratios with stability, since it is derived directly from 
the equations of motion and of heat. Neither theory nor experiment gives 
firm guidance on this point, and it will be assumed that the ratios L,/L,, 
k,, and k,, are nearly independent of stability. The first of these ratios, 
L/L, is simply the ratio of the logarithmic dissipation rates for temperature 
and velocity fluctuations. The large-scale fluctuations of temperature and 
vertical velocity are closely correlated and it seems unlikely that the rates 
of turbulent transfer down the scale of eddy sizes (which determine the 
magnitudes of the conductive and viscous dissipations) could be very 
dissimilar (Ellison (1957) makes an equivalent assumption for a constant- 
stress boundary layer). A more debatable assumption is that the shear 
coefiicient k,, and the temperature—velocity correlation &, are nearly indepen- 
dent of stability, for there are circumstances in which it is certainly not true. 
In the last stages of decay of a turbulent mixing-layer between two streams 
ot different temperature, most of the motion is in the form of gravity waves 


at the interface and permanent mixing on the molecular scale is a rarity. 
3oth these coefficients then approach zero although at different rates. 
Again, in a thick boundary layer of constant stress and constant downward 
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heat transport, buoyancy has a negligible influence near the ground and an 
overwhelming influence at great heights. Here, it is not possible to think 
of k, and k, except as functions of height. ‘The previous analysis has 
assumed stationary, non-developing flow and the use of typical or mean 
values of the velocity and temperature gradients implies that the motion 
in all parts of the flow is essentially similar, so both these sets of conditions 
are excluded. If the flow is everywhere characteristically turbulent, it is 
probable that the correlation factors are always large and near their values 
in constant density flows. 

Returning now to the consideration of equation (3.8), we see that real 
values of the turbulent heat flux and the flux Richardson number are only 
possible if 

, wee, 
. " IZR L, 
assuming, as we do from now, negligible effects of radiative transfer. 
Accepting the substantial constancy of k,, k, and L,/L,, this sets the limit 
to possible values of the flux Richardson number, 
R; < 3. (4.2) 
This limit to the possibility of turbulent motion arises from an impossibility 
of finding any real turbulent intensity to satisfy both the equation for the 
turbulent energy and the equation for the intensity of the temperature 
fluctuations, and it is additional to the limit set by the energy equation 
alone, 


R (4.1) 


R,; < 1. (4.3) 

An interesting consequence of a limiting value for the flux Richardson 
number of one-half is that the turbulent intensity is finite in the critical 
flow. ‘This can be seen by writing the energy equation (3.7) in the form 
(w?)H2 


k, L,/oU ex] 


7 (4.4) 


At the limit R, = }, the turbulent intensity will be about one-quarter of 
the constant density value, and a sudden collapse of the turbulent motion 
occurs as the limit is passed. 


5. COMPARISON WITH EXPERIMENT 
For a proper test of the validity of the analysis, it would be desirable 
to have experimental confirmation that the ratios k,, k,, and L,/L,,, are 
nearly independent of stability, but practically no information is available. 


uw? 


A few measurements of turbulent flows with stable density gradients do 
exist and these may be examined for conformity with the proposed criteria 
for turbulent motion ((4.1) and (4.2)), but nearly all these are of developing 
flows and in these stability is not a simple concept. In an unconstrained, 
developing flow, the Richardson number initially increases with distance 
from the origin of the flow but a comparatively rapid decay of the turbulent 
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motion sets in downstream of the point where the Richardson number 
attains a critical value, this decay corresponding with the sudden collapse 
of the turbulent motion in a non-developing flow. As the decay is not 
instantaneous, the Richardson number may increase beyond the critical 
value. ‘This makes the experimental determination of the critical value 
difficult and measurements of the kind to be discussed can do little more 
than indicate that the critical Richardson number may be near 0-1. 
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Figure 1. Boundary layer on the underside of a heated surface (C. I. H. Nichol). 
Free stream velocity Uy 288 K, 
151 cm, (gx/U%)log(T,,,/T) 1-94. 


150 cm sec |, free stream temperature 7’, 
surface temperature 7’, 388° K, x 


Probably the most accurate and comprehensive measurements in this 
field are those of C. I. H. Nichol in the turbulent boundary layer on the 
underside of a heated surface. In these experiments, an artificially-thickened 
boundary layer was allowed to develop naturally for a distance of 127 cm, 
beyond which point the wall temperature was maintained at approximately 
100°C above the temperature of the free stream. 


At an air-speed of 
150cmsec"! 


, a nearly complete collapse of the turbulent motion was 
observed between 151 cm and 202 cm, the flow then having very low 
values of turbulent intensity and wall stress and a highly inflected form to 
the mean velocity profile (figures 1 and 2). ‘This last effect is a consequence 
of the virtual disappearance of Reynolds shear stresses which earlier had 
kept the air next the wall moving. With their removal, this air slows down 
and the streamlines expand, displacing the outer flow. 


The beginnings 
151 cm (figure 1), although the 
heat trom the wall has not yet spread to the outer parts of the layer, and als: 


of this process are clearly visible at x 
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in the measurements at an air-speed of 190cmsec™! (figure 3). As the 
individual measurements of velocity and temperature have random errors 
reflecting the considerable experimental difficulties, the computation of 
local Richardson numbers is not easy, as is shown by figure 4 in which are 
plotted local Richardson numbers computed by taking finite differences 
between neighbouring experimental points. ‘The scatter is very large and 
only the most docile reader would agree from this that the mean values 
in the outer part of the layers are all less than 0-1. More representative 
numbers have been obtained by using the mean gradients over the outer 
parts of the layers, between 0-5 cm and 2:5 cm from the wall. ‘These 
numbers are given in table 1 and may be more acceptable evidence that the 
Richardson number just before collapse is less than 0-1*. 




















r pr a ee diner 40:04 
Sy ene ceeds Pes 1 y @)" 
Ol | oe P “ oF 
“" | a * ° ; j F mn 
| re 
oe \ 
ae er x 
ma a ' ‘\ 40-2 
U/Up r a | \ 
ra \ log T/T. 
ost/ : . 
a ~ i 40-1 
7, 
L ° 4 See 
7 oe 
ov L j ! ! meen: 
° { 2 ° | 2 3cm 


DISTANCE FROM HEATED SURFACE 


Figure 2. Boundary layer on the underside of a heated surface (C. I. H. Nicholl). 


Free stream velocity U’) — 150 cm sec"', free stream temperature 7) — 288° K, 
surface temperature T°, 388° K, x 202 cm, (gx/U?2)log(T,,/T») 2-60. 


In the second series of experiments, a liquid jet was injected horizontally 
along the interface between a denser and less dense solution of salt in water. 
rom a cine record estimates could be made of the velocity and width of 
the jet. It was found that entrainment of fluid by the jet almost ceased 
when the Richardson number (which increases with distance from the 
nozzle) exceeded ()-3. Beyond this point, the velocity and cross-section 
of the jet remained nearly constant, indicating very little entrainment and 

It should be pointed out that this flow is essentially different from the constant- 
stress layer considered by Ellison (1957). ‘The flow in this developing layer was not 
essentially different from that in one-half of the mixing layer between two streams. 
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Figure 3. Boundary layer on the underside of a heated surface (C. I. H. Nichol). 
Free stream velocity Uy 190 cm sec™, free stream temperature 7, = 289° K, 
surface temperature T',, = 385° K, x = 202 cm, (gx/U?)log(T,,,/To) 157, 
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presumably very little turbulent motion. The first signs of a diminution 
in rate of entrainment appeared at a Richardson number of about 0-05, 
and the critical number, in the sense used in the theory, is believed to lie 
between these limits (see figure 5). 




















Conditions R; Remarks 
Uy 190 cm sec”! = : ae 
en ees mia 0- Collapse just beg ; 
2 = 02cm, T.—T, 96° C. )-070 ollapse just beginning 
U, 150 cm sec-, | 0.046 Spread of heat from the wall not | 


yet complete. 





x tiem, 7 — Po 100° C. | 


Uy 150 cm sec", 0-022 Collapse nearly complete : note 
e = 202 cm, F.— Ty = 100° C. | as) low value of R;. 

















Table 1. 











Figure 5. Spreading of a jet injected along the interface between saline solution of 
densities 1-032 and 1-125 gmcm-*. (a) Variation of widths with distance from 
origin. ‘The straight dashed line is the mean of observations of normal, con- 

bh) Mean velocities of advance of marked fluid put into 


stant-density jets. 
the jet. The dashed line is the mean of observations on a normal jet. 
(c) Instantaneous outlines of the jet, traced from enlargements of ciné records. 
The scale is the same as that of the horizontal scale of the other diagrams. 
(d) Mean Richardson numbers, computed as R; = 0:064g[(p;— p2)/p](D/U;7,), 
where D is the observed width and U,,, is the mean velocity of advance of dye. 


iT 











~~ 


1. Townsend 


a | 
i) 


‘The dithculty in all experiments with developing flows is that a tailure 
of the regulating mechanism will cause the turbulence to decay but, while 
it decays, spreading of the flow and entrainment of ambient fluid continues. 
lor this reason there is a lower critical number at which decay begins and 
an upper critical number attained when decay is complete. ‘The lowe: 
number marks the failure of the normal regulating mechanism of the 
developing flow and should correspond with the one considered in the theory. 

This conclusion that the critical Richardson number is about 0-08 
appears to ignore the frequently quoted series of measurements of the flow 
of fresh-water of salt in the Kattegat, analysed by ‘Taylor (1931) and others. 
In these measurements, Richardson numbers as high as one hundred were 
found although the shear stresses and transport of salinity were much larger 
than could be accounted tor by viscosity and molecular diffusion. ‘The 
observed values of the flux Richardson number were around 0:3. Some 
observations of a similar flow carried out by Dr J. 5S. Vurner in the Cavendish 
Laboratory suggest that this apparent anomaly may be due to the non- 
existence of a critical Richardson number in the ordinary sense for a 
developing flow. In these experiments, fresh water was caused to flow 
over salt water which was coloured for easy identification. Near the 
beginning of the mixing zone, mixing is intense and caused by ordinary 
turbulent movements but, further downstream, the flow settles down to 
nearly laminar flow with a sharp interface between the fresh and the salt 
water. ‘his interface is continuously distorted by irregular gravity waves 
which occasionally break, projecting coloured salt water upwards and, 
presumably, fresh water downwards. Most of these jets seem to fall back, 
losing only a small part of their volume by mixture but clearly transferring 
a considerable part of their momentum by the formation of ‘inviscid’ 
wakes of the vortex sheet type. ‘lransfer of salinity is primarily due to 
molecular diffusion from the jets but is enhanced by the motion of the jets 
through the alien fluid. ‘hese qualitative considerations show that in this 
part of the flow ‘turbulent’ transport of salinity and momentum will each 
be considerably greater than the molecular rates, but that proportionately 
momentum transfer is very much more intense than salinity transport. 
It is probable that the Kattegat measurements were carried out under 
similar conditions and that the measured Richardson numbers do_ not 
refer to turbulent motion but to this random wave motion of the interface. 


[ am indebted to Mr C. I. H. Nichol for permission to abstract from his 
unpublished work the measurements that are represented in figures 1-4. 
| have benefited from observation and discussion of the experiments of 
Dr J. S. Turner on the flow of fresh water over salt water. 
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SUMMARY 


In the case of turbulent flow in a pipe there is a lower 
experimental number to the Reynolds limit for which fully 
developed turbulent flow occurs. Krom the similarity and close 
agreement of the curves showing the coefficient of skin friction c, 
as a function of the Reynolds number R, (based on the momentum 
thickness @) for the circular pipe and flat plate, it is suggested that 
there should be a lower limit to R, for fully developed turbulent 
How on a flat plate. Rather limited experimental data confirm 
this and place the lower limit at R 320. ‘Lhe choice and size 
of transition device is examined in relation to this minimum R, 
and an approximate theory leads to a ‘wire’ Reynolds number 


in fair agreement with experience, 


1. INTRODUCTION 


In the literature it is common to find a curve depicting the turbulent 
skin friction coefficient c, of a flat plate as a function of the Reynolds 
number R., based on the distance x from the leading edge. It is assumed 
that transition occurs at the leading edge and that the boundary layer has 
zero thickness there, i.e. R, = 0. ‘This is an assumption that has a certain 
analytical convenience, but it has dangers for the experimenter and the 
engineer. ‘lhe former is tempted to induce transition at the leading edge 
by means of a wire or other device without giving much thought as to what 
are the real conditions just downstream of it, whilst the latter makes 
estimates of drag for this assumed state which mav be considerably in 
error at low values of R,. 

In this connection the recent paper by Dutton (1956) on the turbulent 
boundary layer of a flat plate brings out clearly the marked effect of the size 
and type of transition device on the distribution of c, and R, considered 
as functions of R,. However, Dutton finds c, to be an approximately 
unique function of R, and independent of the transition device used. 
This is consistent with the assumption of universal ‘inner’ and ‘outer’ 
laws for the’ velocity distribution which has been the basis of important 
papers by Landweber (1953) and by Coles (1954) and will also be used 


in this paper, 
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In this note, with the aid of pipe flow results and the results obtained by (! 
Dutton and earlier investigators on a flat plate, it will be suggested that t 
there is a lower experimental limit to the value of R, for the turbulent ’ 
boundary layer of a flat plate as for turbulent flow in a pipe. ‘Theoretical 
support for this suggestion is given and the functions of a transition device in 
relation to this idea of a lower limit are then considered. 
. 
2. FLOW IN CIRCULAR PIPES i 
Figure 1 shows Nikuradse’s (1932, 1933) results for 7,,/($pU,,) as a t 


function of R= U,,d/v and the roughness ratio d/2e, where U. 
mean velocity in the pipe, d is the pipe diameter, 7,, is the skin friction, 
As R is decreased, the various curves tend 
to that for the smooth pipe (turbulent flow) with the exception of the one 
At a value of log,)R of about 3-65, 


and « is the roughness height. 


for the largest roughness d/2e = 15. 


all the curves, which have now merged into one, travel down a common 
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Figure 1. Resistance of rough pipes (after Nikuradse). 


transition curve which joins the curve for laminar flow at log, R 


This appears to be the maximum value of R for laminar flow with very 
disturbed entry conditions and, even with comparatively large closely 
packed roughness, the skin friction has the laminar flow value at Reynolds 


numbers below this critical value. 
Now the significant point from the standpoint of the present note is 
that it is impossible to obtain results which would extend the smooth pipe 


aoe. 
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(turbulent flow) curve to values of log,,R < 3-65 or R < 4:47 x 108, and 
this is the lowest Reynolds number for fully developed turbulent flow in 


a pipe. 


Defining the momentum thickness for axisymmetric flow by 


paz u 2r 
7 cod | a (1 _ r)( 7) dr, 
o U; U; d 
where U, is the velocity of the pipe centre and wu is the velocity at radius r, 


it is possible to obtain 26/d as a function of the Reynolds number R when 
the velocity distribution is known. ‘This has been done by Ross (1952) 
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Figure 2. Skin friction coefficient against momentum thickness Reynolds number. 


in a re-analysis of Nikuradse’s results to yield R, = U,6/v as a function 
of R, and c, = 7,,/(4pU?) as a function of R,, for turbulent flow. The result 
has been expressed in the form 
cp? = 3-8 logy) Rat 4-4. (1) 
This is shown in figure 2 together with the laminar flow relation 
Ce = 4Ry. (2) 
The minimum Reynolds number for which equation (1) is applicable 


is R = 4-47 x 10°, and the corresponding value of R, is 193. ‘The maximum 
value of R, for laminar flow with very disturbed entry conditions is 87:5, 
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3. "THE FLOW ALONG A FLAT PLATE 

Of a large number of papers dealing with the problem it is proposed 
in this paper to refer to only three, the paper by Dutton (1956) already 
referred to, and the papers by Landweber (1953) and Coles (1954), both 
of which review much of the existing work on the subject. In the 
calculation of the growth of the turbulent boundary and the distribution 
of skin friction, both Landweber and Coles start from assumed velocity 
distributions through the boundary layer. 

In a finite region next to the surface, the ‘inner’ law is assumed to 
apply and the velocity wv is given by 


(3) 


where u_ = (7,,/p)'* and y is the distance from the surface. In the outer 
region of the boundary layer, the ‘ outer’ law applies in the form 
U-u 4 
ui $(/9), (4) 





where U is the velocity at the edge y = 6 of the boundary layer. 

‘The assumption (supported by experiment) that these laws overlap over 
a finite region leads to the well-known logarithmic forms of the equations (3) 
and (4) in this region 


uju, = A+ Blog(u_y/v) ( 


‘st 
— 


(U—u)/u, = C—Blog(y/s). (6) 


‘he ‘inner’ law equation (5) is assumed to be universal for all flows, but 
the ‘outer’ law for the plate differs from that for the pipe in the value of 
the constant C in equation (6). ‘The assumption of the ‘inner’ and ‘ outer’ 
laws leads in both the plate and pipe problems to unique relations between 
c, and R,,. 


In the case of the flat plate flow, the momentum equation 
be, = d0/dx (7) 


enables the distribution of c, and the growth of @ to be found for given 
starting conditions of the turbulent flow. ‘lhe relations between c,; and R, 
as obtained by Landweber and Coles differ because of the different values 
of the constants A, B, C assumed in the analysis depending on the 
individual assessment of the available experimental results. ‘he relation 
obtained by Coles is utilized in this note because Dutton’s results are close 
to it, and in figure 2 it is given for comparison with that for turbulent flow 
in a pipe, equation (1). The relation for laminar flow over a flat plate 


cy = 0-441, ! (3) 
is also plotted. 


The general closeness and similarity of the results for pipe and plate 
is rather striking and it seems reasonable to expect that, as in the case of 
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the pipe, there is a minimum value of R, for turbulent flow along a plate. 
Examination of the results (from various sources) suggests that the minimum 
value of R, for the turbulent boundary layer of a flat plate is about 320 and 
Dutton’s results, which are discussed later, confirm this. 

Now both Landweber and Coles recognized that, in computing the 
boundary layer development along a plate, some initial starting condition 
has to be satisfied as there is aconstant of integration at disposal. Landweber 
in fact restricted his solution to large Reynolds numbers R,, so that the 
starting condition would have negligible effect and he took the constant 
to be zero. Coles fixed the constant by supposing that R, = 0 when R, = 0 
and that the skin friction remains finite. Clearly, the constant of integration 
must be fixed by the value of R, at the transition position and this must 
equal or exceed the value of 320. ‘The value will be fixed by the extent 
of laminar flow and the drag of the transition device. 

Landweber appears to have been the first to recognize that there is a 
limit to the application of the ‘inner’ and ‘ outer’ laws in their logarithmic 
form, since the amount of overlap decreases as the Reynolds number 
decreases. His argument is as follows. ‘The ‘inner’ law equation (3) 
begins to depart from the logarithmic relation equation (5) as the sub-layer 
is approached and the value of y (= y,) when this occurs is often taken as 
given by 

u.y,/v = 30. (9) 


The outer law equation (4) has from experimental data a logarithmic 
behaviour as given by equation (6) up to 

y,/0 : 0-2. (10) 
These equations specify the degree of overlap of the two laws. If y, = y, 
then there is no overlap, and this occurs, from equations (9) and (10), 


when u_6/v = 150 or 
log, )(u_5/v) = 2-176. (11) 


Pipe flow 


We now apply this idea to turbulent flow in a pipe for which Ross (1952) 


C1IVes 
uu, = 5645-6 logy (u,v/v) | (12) 
(1 u)iu ()-785 5-6 log(v 0) 
and, eliminating u, 
Uiu_ = (2/c,)!? = 6°385 + 5-6 log,,.(u_9/v). (13) 


Inserting the value of log,,(u_4/v) from equation (11) into equation (13), 
we find the ‘no overlap condition gives c, = 0-0058, and from figure 2 


log.) R, = 2:305, giving R 202. ‘This limit is marked on figure 2 and 
it is seen that it is very close to the experimental limit for fully developed 
turbulent flow R, = 193. 


F.M. 2B 





378 J. H. Preston 


Plate flow 
Coles’s relations for the velocity distribution in the overlap region are 
uju, = 5°10 + 5-75 logy (u,y/v) | (14) 
(U—u)/u, = 2:80 —5-75 log,(v/5) J’ 
and eliminating u 
Uju, = (2/c,)"? = 7-9 + 5-75 logy (u, 5/v). (15) 
‘The ‘no overlap’ condition gives c, = 0-0048, and from figure 2 the 
corresponding value of R, is given by log, R, = 2°59, ie. R, = 3-89, 
‘The lowest value of R, in the data considered by Coles for which fully 
developed turbulent flow occurs is 320, which is again close to the R, for 
no overlap of the ‘inner’ and ‘ outer’ regions. 

It is not clear whether much significance should be attached to this 
measure of agreement found for both pipe and plate, owing to the rather 
arbitrary method of defining the limits of the logarithmic region by 
equations (9) and (10). Slightly different choice of constants in these 
equations would yield significantly different values of c; and R, for absence 
of overlap. But it is clear that, at values of R, of this order, the sub-layer 
is now an appreciable part of the boundary layer and the viscous stresses 
are becoming important well away from the wall. Experimentally, it 
seems that when the overlap region vanishes further reduction of R, leads 
to the velocity distribution assuming more and more nearly the character- 
istics of laminar flow with a rapid decrease in skin friction to the laminar 
flow value. 

Dutton’s results 


Dutton’s work was carried out at fairly low Reynolds numbers and 
further evidence of a minimum value of R, for turbulent flow over a flat 
plate can be gleaned from his results. 

Figure 3 is a reproduction of figure 7 of Dutton’s report and shows R, 
as a function of R,. for two transition devices, circular wires and sandpaper 
strips. ‘lhe curve for laminar flow has been included in figure 3 and the 
increments to R, caused by wire drag have been estimated by the theory 
given in the next section. ‘The estimated total R, immediately behind 
the wires have been plotted on this figure and are made up as shown in 
table 1. 

In discussing these results we note that R, for the laminar boundary 
layer just ahead of the wire is 120 and that this is less than the value of 
R,, = 162 at which disturbances are first amplified (based on R,* = 420, 
where R,«* is the Reynolds number formed by the displacement thickness). 
It is clear from curve (1) for the 0-013 in. wire that laminar flow is still 
persisting up to R, = 5 x 10° and that fully developed turbulent flow does 
not occur until R,. = 8 x 10°. ‘The repeat set for this wire, curve (2), shows 
fully developed turbulent flow at about R,, = 3 x 10°, transition presumably 
occurring between there and the wire. ‘The difference between these two 
sets has been explained by Dutton as arising from the sensitivity of the flow 
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at the leading edge of the plate to small changes in tunnel stream direction. 
This has probably resulted in the layer approaching the wire being thicker 
than that predicted by the Blasius theory. We note that, on the basis of 
the Blasius theory and the estimated R, for the 0-013 in. wire, R, behind 
the wire is 201 which is slightly greater than that for amplification to occur 
which is R, = 162. Clearly, we cannot regard this value of R, as a lower 
limit to the R, for fully developed turbulent flow, since amplification must 
occur before transition occurs and this does not always occur suddenly. 
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1 @ 0:013° transition wire !" from L.E. (IS set) — L | | 
| j nN — —+ — 13 = 
(2) 0-013" transition wire |" from L.E.(3 set) nie de 
; = results | 
4-0. (3) 0°022. Eransition wire | from L.E = oer Cees! Pe 5 
(4) 2-5" glasspaper strip 3 E 
- 5 2:5 12455D r) 2 
(5) 275 giasspaper strip - 
30 
aXl0 —_+— 
2 











Figure 3. Influence of leading edge conditions on the development of the turbulent 


boundary layer along a flat plate. 














| 

Wire diameter | | Total R, 

a ane | u,/l Udlv AR» | : 
in inches | | behind wire | 
0-013 0-713 423 81 | 201 
——_—_—___ — i ga ee eee real 

| 
0-022 | 0-952 716 | 244 | 365 
| 

Table 1. R, (in front of wire) 120, d = wire diameter, u, value of at y d. 


In the case of the larger wire (diameter 0-022 in.) the corresponding curve (3) 
has fully developed turbulent flow characteristics at the first measuring 
station (R,, = 3:25 x 10°) and it can be extrapolated to the estimated value 
at the wire without difficulty as shown. ‘This suggests that transition occurs 
fairly suddenly at the wire where the estimated R, is 365. Examination of 
Dutton’s curve suggests that a slightly smaller wire than this would still 
have had transition at the wire, so that the previously quoted value of 
R, = 320, as the lowest value of R, for fully developed turbulent flow, 
is in agreement with the above assessment of Dutton’s results. 
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Lhe problem of the minimum Reynolds number for fully developed 
turbulent boundary layers under pressure gradients is more difficult as 
there appears to be no experimental evidence available. If we take the 
agreement between the values of R, for ‘no overlap’ of the inner and outer 
laws and the minimum values of R, observed for fully developed turbulent 
flow in the case of the pipe and plate and being true generally, then it is 
possible to make some predictions. It is known that the outer limit of the 
logarithmic region in terms of the boundary layer thickness tends to become 
less in adverse pressure gradients, so that y,/6 in equation (10) assumes 
smaller values than 0-2. Hence in equation (11), the value of u,6/v for no 
logarithmic region to exist will increase in adverse pressure gradients. 
It is also known that 6/5 increases progressively towards separation, as 
does U/u,, and hence the value of U@/v for ‘no overlap’ of the inner and 
outer laws willincrease. Hence on the basis of the pipe and plate comparisons 
it is expected that the minimum value of R, for a fully developed turbulent 
boundary layer will increase in an adverse pressure gradient and decrease 
in’a favourable gradient. 

We have noted, from figure 3 of this paper or figure 7 of Dutton’s paper, 
the very considerable influence of the type and magnitude of the transition 
device on R, in the range of 0 < R, < 2x 108 and the same is true of 
the skin friction given by figure 8 of Dutton’s paper. ‘This is of importance, 
not only in fundamental experiments which are often carried out at low to 
moderate Reynolds numbers, but also in experiments on ships in towing 
tanks when the Reynolds numbers are of the order of 10‘—2 x 10%. It is 
also important for compressor blading of gas turbines at high altitudes, 
when the Reynolds number may be as low as 10°, and also in experiments 
in small high speed wind tunnels and low speed atmospheric tunnels, 
where Reynolds numbers of the order of 106 are usual. Another point 
of importance is that, even when transition occurs at the wire, increase 
in its size naturally leads to an mcrease in drag and of momentum thickness 
downstream, but this increase in drag is not equal to the extra drag of the 
wire. In fact it will be less because, owing to the increased momentum 
thickness behind the wire, the skin friction will be reduced downstream. 
Thus allowance for the drag of the wire in profile drag measurements, 
for example, by subtracting an estimated wire drag is not correct. 


4. "TRANSITION DEVICES 
On first thoughts it might appear that the primary function of such 
devices is to produce disturbances which will cause transition from laminar 


to turbulent flow. On this basis we could never hope to establish turbulent 
flow nearer to the leading edge than that position at which disturbances 
in the laminar boundary layer become amplified. In the case of the flat 
plate this distance is given by R, = 6x 10', corresponding to R, = 162. 
Clearly any transition device has a drag and as a consequence it increases 
the momentum thickness of the boundary layer flowing over it. ‘The study 
of Dutton’s results in the previous section shows that the drag aspect is 
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all-important when transition close to the leading edge is desired. It 
appears from Dutton’s results that when R, — 320 (the minimum value 
for a completely turbulent boundary layer) then transition is sudden and 
close to the wire. If R, is less than this behind the wire, laminar flow may 
persist to much higher values of R, and transition take place slowly. 

On the other hand, provided R, — 162, thena transition device producing 
vigorous disturbances in the right range of frequencies ought to start the 
transition almost immediately and this should be complete when R, = 320. 
Also for laminar layers with large R, (320), the production of 
vigorous disturbances is the important property. ‘Thus any transition 
device must be considered in terms of both its drag producing and 
disturbance producing potentialities. 

I.et us now consider the drag aspect. ‘he drag ) per unit span can be 
related to the increase in momentum thickness A@ associated with it by 


D = pU?A0. (16) 

In this expression we ignore any contribution by the friction in the immediate 

wake of the device. Now in the case of circular cylinders immersed in thick 

turbulent boundary layers, Sacks (1956) in an unpublished report found that 

Cy = Di(kpuzd) = 0-75 (17) 

over a fairly wide range of u,d/v. In this relation w, is the velocity at vy = d. 
Hence from equations (16) and (17), 

A@/d = iC, (u,/U)*, 


or 


AR, = §Cp(ua/U)* Ry. (18) 


If (R,), is the thickness of the laminar layer in front of the wire and 
if the minimum R, for a fully developed turbulent layer is 320, then 
(R,), + Ry = 320, or 

R, = 320—(R,),. (19) 

From equations (18) and (19), the Reynolds number based on the wire 
diameter is found to be 

R, = (2/Cp)(U/ug)?{ 320 — (Ry); 5. (20) 

‘This equation shows how the size of the transition device is determined 
from momentum considerations. When (R,), is small it is clearly necessary 
to provide this increment AR, in momentum thickness by means of the drag, 
When (R,), ~ 320, then no additional increase in momentum thickness 
is required, but a transition device is still necessary in order to provide 
disturbances which will cause a rapid transition. 

For a wire placed very close to the leading edge so that (R,), = 0, 
(u,/U) = 10 and, taking C,= j, R, U,lv = 3x 1x3 = S50. 
For a wire placed as in Dutton’s experiments so that (R,), 
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we have from laminar boundary layer theory that u,,/U = 0-916, and hence 
R, = 3 x (0-916)? x 200 = 636. A recommended figure, based on 
experimental evidence for producing transition at a wire which is placed 
close to the leading edge, is 

R,, = 600, (21) 


which is of the same order. ‘his is further support for the idea of a minimum 
R, for fully developed turbulent flow. 

When the initial (R,), > 320, the function of the transition device is 
to produce disturbances which will bring about transition as close to the 
wire itself as possible and apart from the rough rule equation (21), it is 
not possible to give further guidance in the selection of the size of the wire. 

Circular wires and other bluff obstacles may, as already mentioned, 
be expected to have a large drag coefficient and to produce considerable 
disturbances in the form of cast-off eddies. As regards the drag, Sacks 
(1956) has found the measured drag of circular cylinders resting on the 
surface in a boundary layer to be given by C,) = D/(Spu5d) = 0-75, which is 
about § of that found for cylinders ina uniform stream. Recent experiments 
by Arie & Rouse (1956) on normal plates with a splitter plate placed 
symmetrically behind give C, = 1:38—1-4, which again is about 3 that 
of a normal plate in a uniform stream. Incidentally, there is in this paper 
confirmation of the relation between D and A@ (equation (16)). On the 
score of drag, the normal plate is to be preferred to the circular cylinder, 
since, for a given AR,, the size of protuberance will be about half that for 
the circular wire. 

As regards the disturbances (when the obstacle is resting on a wall) 
these do not occur as a regular shedding of eddies in the form of a ‘ Karman 
street’. A separation ‘bubble’, some 6-8 diameters in length, is formed 
behind the obstacle. ‘he bubble is composed, in the main, of one large 
eddy in slow motion, but there is a very disturbed region at the tail of the 
eddy, which presumably feeds disturbances into the boundary layer. 

In the case of the circular wire, when the value of R, just behind the 
Wire is just greater than 162 (the value at which disturbances become 
amplified), it seems that the disturbances are either of the wrong frequency 
or are not sufhciently large to bring about immediate transition. In Dutton’s 
paper there is a good example of this where the laminar layer proceeds a 
considerable distance aft of the wire before transition slowly takes place. 


5. CHOICE OF TRANSITION DEVICES 
Circular wires 
Circular wires have found favour because of their convenience. ‘They 
are available in a wide range of sizes and are accurately manufactured. 
‘They are fairly easily stuck to the surface, but care is necessary to ensure 
a uniform spanwise height. From what has been said previously their 
size is governed both by their ability to thicken the boundary layer and to 


shed sufficiently strong disturbances. 
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Normal plate or strip 

‘The use of a thin narrow strip has certain attractions aerodynamically. 
For a given height, it has roughly twice the drag of the circular wire and 
presumably the disturbances it produces will also be stronger. ‘Thus for 
a given increment of AR, it need be only half the height of the wire. ‘The 
fixing would give more trouble than for the wire, as would ensuring 
spanwise uniforraity of height. It might be held in a slot in the surface 


and accurately milled down to size. 


Slotted or toothed strip 

With this device the drag might possibly be raised above that of the 
unslotted strip and the teeth would produce strong three-dimensional 
eddies in the outer part of the boundary layer. ‘This, for fundamental 
work, may be a disadvantage, as turbulence produced in this region might 
influence the form of the ‘outer’ velocity distribution in much the same 
way as excessive stream turbulence. In the practical problems of suppressing 
laminar separation by early transition and of delaying turbulent separation 
by vigorous mixing, they have a definite advantage. 


Sandpaper strip 

This has been used to produce a rapid thickening of the boundary 
layer and as a transition device. It has disadvantages for fundamental 
work in that the results with different strips of the same nominal size are 
not always capable of being repeated and it is difficult to ensure spanwise 
uniformity. Also Dutton (1955) has found that the outer velocity distribution 
does not settle to the same distribution as that which occurs when a wire 
is used. Presumably, this is because eddies are shed from the tips of the 
largest spikes and these travel down the outer part of the boundary layer. 
Air jets 

‘Vhese have been used to promote transition but they cannot be regarded 
as a simple device. ‘here are difficulties in ensuring spanwise uniformity 
and dangers from over-stimulating the outer part of the boundary layer in 
fundamental work. Distributed blowing through a strip of uniformly 
porous material could be effective in thickening a boundary layer. 


6. CONCLUSIONS 
In the case of pipe flow there is a lower experimental limit to the 
Reynolds number for fully developed turbulent flow. From the similarity 
and close agreement between the curves of c,; vs R, for both fiat plate 
and circular pipe, it was anticipated that there would be a lower limit to 
the Reynolds number R, for fully developed turbulent boundary layer 
How on a plate. Rather limited experimental information confirms this 
and places the lower limit at R, = 320. 
‘There is also a lower limit to the Reynolds number for which there is 
an overlap of the ‘inner’ and ‘outer’ laws, the logarithmic region. 
Estimates of the Reynolds number agree fairly closely with the limiting 
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Reynolds tor turbulent flow in the case of the pipe and flat plate. If this 
agreement is to be expected generally, then, on this basis of ‘no overlap’, 
it would seem that the lower limit of &, for fully developed turbulent floy 
should decrease in favourable pressure gradients and increase in adverse 
gradients. 

This concept of a minimum Reynolds number for turbulent flow is 
important when considering the drag of flat plates at medium Reynolds 
numbers and in considering the size of transition device to be used. For 
quick transition near the leading edge, the drag of the device and its ability 
to produce disturbances are both important. An approximate theory has 
been developed on the basis of a minimum &, for turbulent flow which, 
in the case of circular wires, leads to a ‘wire’ Reynolds number in fair 
agreement with that suggested by experience. 

As regards the choice of transition devices, the conclusion is that a 
normal strip has the least height to achieve transition, but the circular wire 
is more convenient. When suppression of laminar separation and delay 
of turbulent separation is important, it is thought that a toothed spoiler 
or slotted strip should give good results; but the generation of eddies in 
the outer part of the boundary layer is objectionable in fundamental work, 
as the ‘outer’ part of the velocity distribution may be affected. 
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SUMMARY 


Investigations are made of the plane flow of an inviscid 
isentropic gas at high subsonic and sonic speeds past a finite 
wedge of small angle set at zero incidence in a channel with parallel 
walls. Hodograph methods are applied to the determinatidn of 
the stream function, which, under the usual transonic approxi- 
mation, is a solution of ‘Vricomi’s equation. In § 2a brief summary 
is given of the derivation of the fundamental solution of this 
equation in terms of Bessel functions. 

‘l'wo models of the flow pattern are discussed. ‘lhe model of 
Cole is examined in $3 with the sonic line from the shoulder 
extending completely across the channel at right angles to the wall. 
In §4 the Helliwell—-Mackie model is taken in which there is a free 
stream breakaway at sonic velocity from the shoulder of the wedge 
and the velocity far downstream may be either uniformly subsonic 
or sonic. Values are obtained for the drag coefficient in both cases 
and a high degree of both qualitative and quantitative agreement 
is found. On the basis of the free streamline theory explicit 
formulae are given relating, in terms of the channel width and 
upstream Mach number, the drag of the wedge in the channel to 
that of the same wedge in the free stream. 

Attention is drawn in §5 to certain properties of plane jet 
Hows that may be deduced from the investigation of the flow 
with free stream breakaway past a wedge in a channel. 


1. INTRODUCTION 
‘The two-dimensional steady flow of an inviscid isentropic gas at high 

subsonic and sonic speeds past a finite wedge profile in a free stream has 
been investigated by many workers, employing various methods based upon 
the hodograph transformation by which the resulting differential equation 
for the stream function is made linear. A common procedure, which 
will be followed in the present paper, is to use the so-called transonic 
approximation in which this differential equation becomes the well-known 
‘Tricomi equation 

ae oy 

- u=~—> 0, 


A > 
cu” Cv" 

































386 ¥. B. Helliwell 


where U = a*(1—u/(y + 1)) and V = a*v/(y + 1) are the cartesian components 
of velocity, the axes being taken parallel and perpendicular to the uniform 
flow at infinity, a* is the velocity of the gas when the Mach number is 1 
and y is the adiabatic index of the gas. In this approximation V = () at 
infinity upstream of the wedge and u, v are small. Henceforth the suffix 1 
will refer to conditions at infinity upstream. Furthermore, it is known 
that ‘’ is proportional to the y-coordinate. 

It should be remarked* that recent work by Vincenti, Wagoner & Fisher 
(1956) shows that even more accuracy may be expected from the transonic 
approximation if solutions of ‘T'ricomi’s equation 
224" q2y” 
=" 








a +t 

ou? dv? 
are used in which uw ts interpreted in a somewhat different manner from 
that in the present paper, and ‘” is the stream function ‘’ multiplied by 
a certain other function of w. 

Guderley & Yoshihara (1950), Cole (1951) and Helliwell & Mackie 
(1957) have obtained solutions of ‘Tricomi’s equation under different 
hypotheses regarding the flow pattern downstream of the shoulder of the 
wedge. In each model the flow, accelerating from zero velocity at the tip 
of the wedge, attains sonic velocity at the shoulder. ‘The solution of 

Guderley & Yoshihara is of considerable complexity since it includes a 
determination of the flow field between the sonic line and the limiting 
characteristic of the Prandtl-Meyer expansion at the shoulder. Less 
difficult analysis is involved in Cole’s solution in which the supersonic 
region degenerates into a sonic line straight and normal to the direction 
of the flow far upstream; the flow pattern is investigated upstream of this 
line. In the work of Helliwell & Mackie the sonic line is specified as a 
free streamline starting from the shoulder and the velocity is subsonic or 
sonic throughout the entire field of flow. A comparison of these various 
models shows a close similarity in the pressure distribution over the wedge 
nose upstream of the shoulder and hence correspondingly small variations 
in the estimate of the drag coefficient of the wedge. 

An extension of the approach of Guderley & Yoshihara to the situation 
when the wedge is placed symmetrically in a wind tunnel has been carried 
out recently by Marschner (1956). However, he only investigates the case 
of the choked tunnel. In the present paper similar investigations are made 
of the flow past a wedge at zero incidence in a two-dimensional channel 
with parallel walls under more general conditions of downstream flow. 
In §3 the basic model of Cole is used with the sonic line extending from 
the shoulder of the wedge normally to the channel wall. In §4 the 
Helliwell-Mackie model is taken with the sonic line appearing as a free 
streamline from the shoulder. Expressions are obtained for the drag 
coefficient in each case and the channel walls are found to produce similar 
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etfects upon its value, as compared with that in a free stream, for both models 
of the flow. In the free streamline model, equation (25) taken in con- 
junction with figure 8 gives correction terms which yield explicitly, in terms 
of the upstream Mach number and the channel w idth, the relation between 
the drag of the wedge in the channel and of the same wedge in a free stream 
with the same Mach number. When using these results it should be ensured 
that the channel width has been taken sufficiently large for the flow far 
downstream not to become supersonic; the limiting relationship between 
channel width and upstream Mach number for uniform sonic velocity 
downstream is shown in figure 6. 

Finally in $5 attention is drawn to a result concerning jet flows that 
may be deduced from the investigations of the flows with free streamlines 
past the wedge in a channel. 


2. FLOW PATTERNS IN A SUBSONIC STREAM 

In the subsequent analysis we shall be concerned with flow nowhere 
supersonic. In this section we summarize briefly the equations relating 
the variables which define the flow patterns. ‘These relationships only 
hold, of course, within the limits of the transonic approximation. 

‘he dimensionless velocity variables u and v have already been defined. 
The Mach number M is given by 1—M? =u. ‘The polar angle 6 of the 
velocity is related to v by v = (y+1)@. ‘The local pressure coefficient is 
defined by 

- M1, 172 
, = (P—p,)/(se1 U7), 
where p is the pressure and p is the density. It is straightforward to show 
that, according to the linearized transonic theory, 
6, = Au—m) [ly +1). (1) 
‘The equation of continuity and the condition for irrotational flow are 
respectively 


uu Tad UV _—— Q, ,, £. 2, — 0, 


which may be inverted to yield the hodograph equations 


If x is eliminated these lead to the equation of ‘Tricomi 
Vuut UV yy = 9, 

in which u > 0 represents subsonic flow. If we set 

r= 343? = 2(1 — M?)32, (3) 
then, on solving ‘T'ricomi’s equation by the method of separation of variables, 
simple solutions are obtained of the type 

y = Ae@ Or), (4) 
where @ , ,3(Ar) is any linear combination of Bessel functions of the order 
indicated and A is any constant, real or imaginary. 
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3. MopeEL |. SONIC LINE NORMAL TO CHANNEL WALL 


In this section we employ the model of Cole for flow at near sonic 
velocity past a wedge profile. ‘The basic feature of this model is that the 
sonic line from the shoulder is straight and normal to the flow at infinity 
upstream. It is pointed out by Cole in his original paper that, so far as the 
flow over the front portion of the wedge is concerned, it is a good approxi- 
mation to replace the actual shape of the sonic line in this way, and hence 
a good estimate should be obtained for the drag coefficient. At low subsonic 
Mach numbers the supersonic region is small and so the approximation 
is close to being correct; for larger Mach numbers, whilst the actual sonic 
line is by no means straight away from the wedge, yet the upstream influence 
of the errors in the supersonic region introduced by the approximation 1s 
weaker. 
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Figure 1. Model 1. Physical plane. 


Consider then a wedge of semi-angle 6 placed symmetrically in a channel 
of semi-width A as shown in figure 1. The x-axis is taken as the axis of 
symmetry and the y-axis passes through the tip of the wedge. Hence, 
we need only consider the flow in the upper half-plane y > 0. The channel 
is supposed to be blocked by the sonic line which extends normally from the 
channel wall to the shoulder of the wedge. This is represented by BC. 
AOB is the dividing streamline ‘? = 0 on which at A, the point at infinity 
upstream, the velocity is subsonic with Mach number M,, at O, the wedge 
tip, the velocity is zero and along the wedge face OB the flow accelerates 
to attain sonic velocity at B, the shoulder of the wedge. DC is the stream- 
line 4 =k with associated vy = K corresponding to the channel wall along 
which the flow accelerates steadily from its subsonic upstream velocity 


to sonic speed at C. Along the wedge face OB we also have v = vy = (y + 1)6. 
‘The boundary value problem may now be established in the hodograph 
plane, as shown in figure 2. CD is the line y = K. AOB is the line y = 0. 


PBC is the sonic line x = constant, from which it follows, since u = 0, that 
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x, = 0, whence using (2) we have y, = 0. ‘I'wo other lines of constani 
are shown. Clearly there is a finite singularity at the point DA. 

‘The problem is now stated in terms of rand v. ‘The condition along BC 
becomes r!%y,, = 0 when r = 0 for all 0 < v7 < a which is automatically 
satisfied provided that y, remains bounded atr = 0. ‘The simple solutions (4) 


show that we must take y, = Oatr = Oas the condition to fit this requirement. 





Hence we have 
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Figure 2. Model 1. Hodograph plane. 


‘The linearization principle, applied to the stagnation condition at the tip O, 
leads to 
c= 0, y = 0, as 7 —> oO. 


Finally, taking the wedge to be of unit length, 


r=) atv = % r= 0. 








A solution of the problem is obtained by taking a combination of the 
simple solutions (+) in the form 
-. : sinh A(vy — 7) 
E, | (Ayr? 8J_, (Ar) - ~——" d), 


Pe sinh Az, 


wn 
— 


Differentiating under the sign of integration, we have 





oo z sinh A(vy — v 
V — | XA f(A)r?3.J5(Ar) —= (o 
| : / 0 
Here f(A) is an arbitrary function of A that is to be determined so that (5) 
satisfies the first two boundary conditions, since it is immediately apparent 
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that the last two are satisfied automatically. Hence f(A) is given by 


- Oo Cl 
| f(A)r'3J_,(Ar) dA = K, c= < Fy, 
|) r>1y. \ 
An application of the Hankel inversion theorem and a simple integration 
then yields the result that 
fA) =A KPBS (Ar) dr = Kr? J y5(Ar). 
“0 | 
Inserting this in (5), the solution of the problem becomes ¢ 
sinh A(v» — v) 
— Kp3yil3 TOP eee et 
y = Kr | I yg(M4) 113A") Xe, dn. (6) 


‘That this form has a singularity of the correct type at r = r, may be confirmed 
in the usual way. It may also be noted that (6) simplifies into Cole’s solution 
for the wedge in the sonic free stream if Kr}’* tends to a finite constant as 
r, > 0 so that the channel becomes infinitely wide as the upstream Mach 
number tends to unity. ‘his behaviour of K is confirmed in a later result 
(equation (9)). 

The x-coordinate is now obtained from equation (2). After a little 
analysis we find that 


cosh A(v») — v) 


sinh Avy 


ee L(g Ke | J 33(Ar)1 23(Ar) 


~ oO 


dX. (7) 


‘The unit additive constant in this expression arises from the condition that 
the sonic line on which r = 0 is x = 1 when the wedge is taken to have unit 
length. 

In the subsequent work a series representation for the x-coordinate on 
the face of the wedge will be useful. ‘lo obtain this we set v = vp in (7) 
and expand cosech Avy in partial fractions by 


oO a 1 "ry 
cosechAvy = — +2 — 


Av n=1 (Avo)? - (nz)? ; 


After some manipulation and use of standard results concerning integrals 
involving Bessel functions, we ultimately derive the series 
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Because of the asymptotic behaviour of K,,(s) for large z, the stagnation 
condition, x = 0 as r -» %, applied to the latter representation gives 
7 A/3 2\4/3¢, 
Kr}? = 2(%)*%u,. (9) 
Writing 7, in terms of M, through (3) and replacing vy, by (y + 1)6, we thus 
have the result. 
[3(y + 1) 8K = 2[(1 — M3){3(y + 1)}24]. 
‘This is the relation between the semi-channel width and the upstream 
Mach number that, on account of continuity of the flow, holds for the 
present model, if the channel is of such a width that sonic velocity is attained 
everywhere across a section through the wedge shoulder. 
‘The drag coefficient is defined in the customary way as 
Cp = Di(2p, U7), 
where D is the drag on the upper face of the wedge and, according to the 
linearized theory, is given by 
D = 8 C, dx, 
the integrations being taken over the length of the wedge. Introducing the 
relations (1) and (3) it can then be shown that 


B)23§ ( ~~ Ay 
Cp = 2) T 1 | (=) dr+r 
y+1 [Jo ae 
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(10) 


The series representations (8) may now be substituted into (10) and after 
considerable algebra, in which (9) is used to replace A in terms of r,, an 
expression for C, in the form of an infinite series is obtained as 


Be! Cy _ a (= ” s < csr K; (= i 
om ka (§)x' e 2v9 m2 al o vy 





where ¢(z) is the zeta function of Riemann. For values of the upstream 
Mach number MM, little different from unity, a series expansion in powers 
of (1— M7) may be obtained for C,, in which the leading terms are 

26 ‘ _ 

— (1- Mj) +001 — iy, (11) 
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Y 1* 
Cp = Cp- 


where C7, is the finite value of the drag coefficient in Cole’s model for the 
wedge in the sonic free stream. As M, decreases from unity, the value 
of Cy for a wedge in a channel diverges from the corresponding value of C,, 
for an identical wedge in a free stream with the same upstream Mach number, 
but reference to equation (72) of Cole’s paper shows that the difference is 
of the order (1—M?)*. A specific comparison is shown graphically in 


figure 3, 
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+. Mover 2. SONIC LINE A SECTION OF A FREE STREAMLINE 

We now turn to investigate the flow pattern past a wedge profile symmetri- 
cally placed ina channel, employing the model incorporating a free streamline 
proposed by Helliwell & Mackie. For the wedge in a channel the following 
general properties of the flow are laid down. ‘The dividing streamline ‘VY = (0) 
is straight from infinity upstream, where the velocity is subsonic with 
Mach number .V/,, to the tip of the wedge where there is a stagnation point. 
‘The upper branch of 4" = 0 then follows the upper face of the wedge, the 
gas accelerating to attain sonic velocity at the shoulder at which point the 
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Figure 3. Miodel 1. Drag coefficient. 


streamline breaks away from the wedge, retaining sonic velocity until it 
again becomes parallel to the channel walls. ‘hereafter it remains straight, 
the velocity of the gas along it meantime decelerating from sonic speed 
to some high subsonic velocity with Mach number WM, (> M,). ‘he 
critical case of a channel with uniform sonic velocity everywhere in the flow 
field far downstream is merely a special case of the above with M, = 1. 
\lore generally, the model describes a tlow pattern with uniform subsonic 
velocity upstream, a local sonic region in the neighbourhood of the shoulder 
of the wedge and ultimately a uniform subsonic velocity downstream in 
addition to the wake, regarding which certain information may be deduced. 
Irom a theoretical point of view there is the added advantage that throughout 
the whole flow field the velocity is nowhere supersonic and hence limit lines 
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cannot occur to render invalid the transformation back to the physical 
plane. 

‘The pattern described in the preceding paragraph is shown in figure 4. 
The semi-channel and semi-wake widths are K and H respectively. ‘lhe 


upstream portion of the dividing streamline ‘’ = 0 is taken as the line of 
the x-axis with origin at the tip of the wedge. ‘The channel walls are thus 
the straight lines y + K. Because of symmetry only the flow in the upper 
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Figure 4. Model 2. Physical plane. 
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Figure 5. Model 2. Hodograph plane. 


half-plane y > 0 need be considered. AOBCD, then, is the dividing 
streamline. At the point A far upstream the Mach number is M, with 
associated r = r,; at O the velocity is zero. Along OB, the wedge face, 
the component,v = vy) = (y+ 1)8, where 6 is the semi-angle of the wedge; 
the gas accelerates to sonic velocity at B and retains this speed along the 
section BC. At C the velocity is again parallel to the x-axis and along CD 
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the gas decelerates until at D far downstream the Mach number is M, with 
FE is the upper channel wall along which the 
gas velocity varies from Mach number ™, at F to M, at E. 

The boundary value problem is now set up in the hodograph plane which 
is shown in figure 5. It is recalled that, to the order of the transonic 
approximation, ‘’ is proportional to vy. AOBCD is the line y = 0, EF is 
the line y Finite 
singularities occur at the two pairs of coincident points AF and DE, corre- 
sponding to upstream and downstream conditions, respectively. The 
remaining boundary conditions, in terms of r and v, are 


associated r = ry (< 74). 


- K, and two other lines of constant y are shown. 





y= 0 v= U, TE a as 


y= 0, 2) = Q, 7 V5 
y = Q, Z Uo) PSU. 
y = 0, (0) OSU oS 


‘he stagnation condition at the tip of the wedge is, according to the transonic 
approximation, 

6=V, y=U as r-> oO. 
For a wedge of unit length, we also have the further condition 


ee is C=O, ator =U. 





The solution of the present problem may be derived from a combination 


of the simple solutions (4) of the type 


— sinh A(z v) . 
y | g(A)r?3J,,(Ar) — ile, dn. 


Jo sinh Av, 
Whereas this automatically satisfies the last two boundary conditions fot 
any g(A), it only satisfies them all if g(A) is a solution of the equations 


L 


| g(A)r*8Sq)3(Ar) dA = 0, Oar < fy, 


= 0, r>t1}. 


The function g(A) may be determined by an application of the Hankel 
inversion theorem and we obtain 


(A) = K{r3° J _ai9(AY2) — 179° _2/3(A"y)}- 


Hence the formal solution of the problem is 


“2 
> 2 9/2 
y = Kr? | {r5 


J _ajg(AV2) = 17° I _a19(AM)}-F yar) 


| sinh A(vp — v) dv. (12) 


sinh Av, 
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This has singularities at y = 7, andr =r, whenv=0. That they are 
of the correct type may easily be seen by an investigation of the form of y 
in the neighbourhood of r = 7,, r =r. We also note that if we let r, > r, 
and K > o simultaneously, we derive the special case of the flow in a free 
stream. In fact, by carrying through this limiting process on yin equation (12) 
using the relation (14) between K, 7, and 75, we find that we recover the 
solution for flow past a wedge of unit length in a free stream given by 
equation (3) of the paper of Helliwell & Mackie. 

The x-coordinate of the solution is next determined. <A direct sub- 


stitution from (12) into (2) and the resulting integrations yield 
es ’ cosh A(v, — v) 


sinh Av, 


da. 


(13) 


The constant which appears as a result of the integrations may be shown 
to be zero as a consequence of the stagnation condition at the tip of the 
wedge. 

The wedge will be taken to be of unit length and we have therefore x = 1 


at the shoulder where v = vy as r—> ©. ‘Thus, from (13), we find 
P= | 59 (Ag) — 17° J _/9(Ar,)} AA. 


On expanding the Bessel functions, setting Avy = ¢, and interchanging 


orders of integration and summation, it follows that 





($)'3K < (—1)" af V2 2n—4/3 waft \? 1/3 “0 g2n—4/3 Jy 
| = oh » “sri 5 i) { eo eT \ ey ] a Sr er ae 
P(A)e 2! (n+3) * \Sro, 209] sinh ¢ 


[his is simplifte d by the use of a well-known integ1 il representation oft the 


Riemann zeta function ¢(z), giving 


(21 Ky! =e 1 
‘ Al) ay: i} ‘ )2 , 
i SING Ey > Sor | (71 ray te om 1 )C( i) )X 
V7TL (3) n 





This equation determines the relationship between upstream and down- 
stream Mach numbers VM, and M, for a wedge of unit length placed in a 
channel of semi-width K, where 


r, = 3(1—M2y2, ory = 3(1— M2)". (15) 


Figure 6 shows this result graphically in the important and limiting case 
of a channel when the velocity is uniformly sonic downstream (r, = 0). 
For purposes of comparison the corresponding result (9) for the model 
of the previous section is also shown. 

In the subsequent work series representations for the coordinates will 
be required. We shall obtain these by means of intermediate contout 


2C2 
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integral forms. First let us investigate the y coordinate. Setting Avy = t 
the result (12) becomes 


K ay! 3 f a tmro = tor Se 
; 4 V2 *d 2/3 ( *) ry | 9 (=) ; x 
Cy 0 q "y Uo 


KJ 13 a oon dt. (16) 


here are three distinct cases of this integral to be discussed, depending 
upon the value of r. 
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Figure 6. Upstream Mach number and critical channel width. 
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where W, is the part of the integral along the imaginary axis from 0 to1a, 
and set v = fe’"?, vy = te"? in W, and W,, respectively, we find after some 
manipulation that 

gu SW. 


Since the poles of the integrand of W occur at the integer values of v, the 
residue theorem of Cauchy applied to .4W will yield a series form for y. 
We obtain 


2Kris ~ /( nor. nmr, \ ) nr \ nTv 
y = — 2 47K, (“=2) -1°Ke (2) rl, (* )sin(“*), (17) 
ca) 1=1 | . ca) } Vo Up 


In the case when r >7, > 7, considerations similar to the above 





may be used to show that 


iKr' oF ( ans vars ior {vary ae var err —v/ V9) 
ae | {13°93 - 17° sig > Kis — } ———— vv, 
; UM Jol ‘ Up i Uo Up SIN vr 


where C is the same contour as before. ‘The Cauchy residue theorem 








applied to this expression gives 


ae = J 2s NTT » ii nnr,\\,- (nunr\ . (nv 
Ps ee BA a eo oe ee ee 
“CO n=1 { “% i. “6 f Co ar) 


(18) 








Over the range r, <r < r,, the two forms of y already obtained suggest 
that we consider 


- iKri3 ¢ { a VII's \ _, var \ reise vary var \ | 
W' = ——] 4 r3°I_sl — }Ki al — )4 ri’ Kosi — ]L 3 — \ 
Yo Jol ca) cary ca) lo J) 


ela(l—v Ue) 








x ———— dv, 


sin vr 


with C the contour described previously. Once again, the contribution 
from the infinite semi-circle is zero, but we now find that the imaginary 
part of the contribution from the origin indentation is 


— K(1—v/z). 
Hence, proceeding as before, we discover that 
y = SW’'+ K(1—o/v 9), 


and an application of the residue theorem of Cauchy yields 


2Krii® 2 ‘nar nar 
ati } 2/3 ROE ex 7 
a ee z ae | zal = is (= ) t 

“oO "ES at) ba |) 

nar n7r\ (nv v 

° 4 - 1 . - 
-+ rr Kes( ee -sin( =) t K(1 —_ *), (19} 
cai) Uo \ %o \ 9 


‘lhe required series representations for y are given by equations (17), 
(18), (19). The discontinuities when v #0 at r=r,, r=ry are only 
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apparent; continuity may be shown to exist at these points by use of a 
result concerning Bessel functions (Erdclyi 1953) and the Fourier series 
representation 





sin — , v £0. 


Of course, the value of y is still undefined at r =17,, r =r, when v = 0. 

Series representations for x may now be obtained by the use of 
equations (2). ‘The constants of integration are determined by the 
stagnation condition at the tip of the wedge and by the necessary continuity 
of the wx-coordinate everywhere. ‘The three forms corresponding to 
equations (17), (18) and (19) are the following. For the range 0 <r <r, 











Vi» 
2(3)18Kr23 2 .... (nares Sink ite Ne 
7 “( ——— 3 <r ‘Kya 3 2) — 17" Ky3\ — " rx 
ca) L | i 0.7 “0 ; 
= 34/3 
177) NTU >) K 
x | 1s (cos )+ ( = (r5 —15 ); (20) 
Uy carn “Cy 


for 7 >, > 71, 








and for the range rg < r < 1, 


2(2)'8 Kris are nnro\_. (nmr 
Af: = > 772°l 2 (2) Ky { ) 7 
Cp | Up 7 ca) 





We are now in a position to calculate the drag coefficient, Cp, which, 
as indicated in the previous section, is given to the order of the transonic 
approximation by equation (10). ‘The evaluation is quite straightforward. 
The series representations (20), (21) and (22), with v = vo, are inserted into 
the infinite integral and orders of summation and integration are 
interchanged. ‘The resulting integrals involving Bessel functions are all 
standard tvpes and after their evaluation a little algebra leads to the 


exp! eSsion, 


y \2 , )( 3)2/355/3 » \2/3 
ee - '5\* LAS) O (iy G »4 
Cy = 3Ka4 (71) (72) | Aet  aol (23) 
i j 4 
U,/ \Uq/ \Y L)} 


] I F j } 72 1 
For a wedge of unit length, we « xpress A 1n terms or7, andr, by equation | 1+) 
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and thereby obtain 


(y 1 1)'8C,, = ()" : 


«Sb \2ro 








>Y [(-— 1)" tt/nl]3-2" Pn — §)(22"-13 — 1)0(2n — $){(371/209)" — (372/20 9)?" $ 
n=] 
This may be written in terms of the upstream and downstream Mach 
For flow sufficiently near sonic this may be expanded 


numbers M, and M,. 
The leading terms 


as a ‘l'aylor series in powers of (1 — M7) and (1— M3). 
of such an expansion are 
26 " chin pe 
C = €2"- apes (1— M?)+ Of(1— M?)8, (1— M3)3}, (24) 
y+ 


D 


where C},* is the finite value of the drag coefficient of this model for the flow 
past the wedge in a sonic free stream. The result (24) should be compared 


Free stream 
, , Channel! 


+O 0:8 06 04 O02 O 

Figure 7. Model 2. Drag coefficient (critical channel). 

with that for the previous model as given in equation (11). The special 
case of the drag coefficient for a wedge in a critical channel with variable 
subsonic upstream velocity and sonic velocity along the entire free streamline 
and uniformly far downstream is shown in figure 7. ‘This should be compared 
with figure 3 for the corresponding case of the previous section. The effect 
of the channel walls is seen to be very similar for the two models of the flow. 
A more useful form for the drag coefficient in the general case is, however, 
obtained if the asymptotic expansion of equation (23) is taken for large K. 
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The drag coefficient is then written as a series expansion in ascending 
powers of K~! where the relation (14) has been previously used to express r, 
in terms of K and 7,. After considerable algebra we obtain an expansion 
of the type 
| 1,3 IV 
(y +1) iC , aM)  (M,) n c(M,) 





rE cise (25) 


67/8 L as Pehanneld K . KE _ 
where K = K(1—M7)3/[8(y + 1)}>°. The coefficients a, 6, c, ... are extremely 
complicated and inelegant functions of M, which increase slowly as M, 
decreases from unity. The first two are shown graphically in figure 8. 
Cd. ANd Cp... ., denote the drag coefficient of the same wedge in a 
free stream and channel, respectively, with the same upstream gas velocity. 
The expression (25) is thus essentially one yielding correction terms to be 
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Figure 8. Model 2. Drag correction coefficients. 


applied to the drag coefficient in the channel in order to obtain the corre- 
sponding value for the drag in the free stream. It is remarked that the 
correction to be applied is, in general, small within the range of validity 
of the present theory which does not apply unless the channel is wider 
than the length of the wedge by a factor of the order of 10. ‘That this is 
so may be noted from the fact that, for a given value of M,, the value of K 
yielded by figure 6 is that corresponding to uniform sonic velocity across 
the channel far downstream. In general, however, the velocity will be 
subsonic and considerations of continuity of the flow show that the value 
of A thus determined is a lower bound. 
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The location of the end of the sonic line is the point C in figure 4. Here 
the wake becomes parallel to the channel walls. ‘The x-coordinate of this 
point is obtained by setting r = 0 and wv = 0 in (20) and has the value 





2(3)8K 2% (nr\-23 ( ,.,, (nar, since /MBT;\) 
Xo = TWH: —— Zz (=) 12" Kaal es ) ry, Ky,(“4) rt 
(43)% n=1 \““%O \ ©o . %O J 
($)*8K ,, ; 
- - 26) 
2a (r| Yo ) (2 
=o 


This is finite for all 7,, r. 40 so that the wake becomes parallel at a finite 
distance downstream of the wedge, as was laid down in the flow pattern 
for this model. It is of considerable interest to note that this result remains 
true when rv, = 0, so that uniformly sonic conditions in a parallel flow 
would be attained at a finite distance downstream in the case when the 
channel width is critical. 
Finally, the semi-width of the wake can be calculated. From figure 4, 

for a wedge of unit length, we have 

os os 

H-8 =| “ dy = | ° Odx, 

a aoe 

where 4 is here the slope of the sonic line. Integration by parts leads to 


“1 
0 


y+)lH =| xd, 
where w is to be taken along BC. The appropriate series representation 
of x to be inserted in this integral is the limiting form of (20) as r— 0. 
Thus, after the integration is performed, we find 
(3)°K | ' 
H = ~~ (ri3— rh), 
Ay+1)° ' 


which by means of equations (15) may be written 
_ (M2 — M})(2-— M?2- M2 Pp 
Gono it lal (27) 
2(y +1) 
It should be recalled that this value is not independent of 6, as appears at 
first sight, since K, M, and M,, are interrelated with 6 through equations (14) 


and (15). 





5, JET FLOWS 

It is apparent that if we redraw figure 4 with FE as axis of symmetry 
and take the flow region bounded by the dividing streamline AOBCD and 
its image A’O’B'C’'D’ in FE we obtain figure 9. The physical problem 
portrayed is then that of a two-dimensional jet flowing from a convergent 
nozzle through a region of constant pressure into a duct with parallel walls. 
Far upstream in the channel of semi-width A, the Mach number of the 
flow is M,. At the orifice BB’ sonic velocity is attained at the boundary, 
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and the gas issues as a jet with free boundaries BC, B'C’ along which sonic 
velocity is maintained. At CC’ the gas enters a duct with straight, parallel 
walls in which the flow decelerates to a subsonic velocity with Mach number 





M, (> M,). 
bs 
J 
O 7 : we - D 
i B 
x 6 


Figure 9. Physical plane for jet flow. 


Figure 10. Downstream position of parallel sonic jet. 


‘The case when ™, | is of particular interest since then the duct 
may be dispensed with and we have the problem of a two-dimensional jet 


flowing through a convergent nozzle into a region of constant pressure 


which has the value associated with the sonic velocity. ‘The flow in the jet 
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has a uniform sonic velocity far downstream, where the ultimate semi-width 
of the jetis J (= K—H). The results of the previous section in the special 
case when r, = 0 are thus appropriate for the solution of this problem. 
We take FE as x-axis and OO’ as y-axis. Then, as has already been 
remarked, the jet becomes parallel at the point C located a finite distance 
downstream of the orifice. For a nozzle of unit length the value of xo, 
the x-coordinate of C, is given by equation (26) with r,=0. We recall 
that the relationship between the channel semi-width K and upstream 
mach number M, is shown graphically in figure 6. ‘The variation of x 
with M, is shown in figure 10. Finally, it follows from equation (27) that 
the value of J, the semi-width of the jet, is given by 


J (1— M2) 
K-  &Ay+1)" 


As remarked at the conclusion of the previous section, it should be noted 
that this result is not independent of 6 since K and M, are interrelated 
with this angle. 
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SUMMARY 

‘The slow uniform motion, after an impulsive start from relative 
rest, of a paraboloid of revolution along the axis of a rotating fluid 
is investigated by using a perturbation method. ‘The principal 
purpose of the note is to illustrate the mechanism by which the 
Huid is not subjected to any substantial radial displacement, which 
is a direct consequence of the requirement that the circulation 
round material circuits should be constant when the perturbation 
velocities due to the motion of the paraboloid remain small. It 
appears that the mechanism is an oscillatory one in which the 
distance between any fluid particle and the axis of rotation oscillates 
sinusoidally in time with small amplitude. As time progresses, 
the amplitude of the oscillation decays to zero everywhere except 
on the paraboloid. ‘The ultimate motion is then a rigid body 
rotation everywhere except on the paraboloid and the axis of 
rotation, where the perturbation velocities continue to oscillate 
indefinitely with small amplitude. 


1, INTRODUCTION 

‘he motion of bodies in a rotating fluid has been a subject for a series 
of investigations in recent years. ‘The perturbation caused by the motion 
of a body in an inviscid fluid exhibits different characteristics according as 
the fluid is at rest at infinity or is rotating about an axis there. Thus if the 
Huid is at rest at infinity, the flow is everywhere irrotational and dependent 
only on the instantaneous velocity of the body. But if the fluid is rotating 
about an axis, the perturbation in the fluid velocity depends not only on 
the instantaneous velocity of the body but also on its past history and is 
in general neither steady nor irrotational; and even in cases where a steady 
solution of the governing equation can be found, there is no guarantee that 
the flow can be set up by starting the body from rest relative to the rotating 
system. For these reasons it is necessary to consider an initial-value problem 
while dealing with this type of fluid motion. 


When the body moves slowly it has been customary to use a small 
perturbation theory. Using this method Stewartson (1952, 1953) has 
investigated the slow uniform motion, after an impulsive start from relative 
rest, of a sphere and an ellipsoid along the axis of a rotating liquid. In both 
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these cases he found that ultimately the fluid inside the circumscribing 
cylinder C with its generators parallel to the axis of rotation is pushed along 
in front of the body as if it were solid, while outside the cylinder there is a 
shearing motion parallel to the axis of rotation. ‘There is also a swirling 
motion about the axis inside the cylinder C. ‘The ultimate velocity distri- 
bution in the fluid is in general steady and two-dimensional (in the sense 


that the motion is the same in all planes perpendicular to the axis of 
rotation) everywhere except on the body and on its axis, where it oscillates 
finitely. In fact the linearized equations show that every slow and steady 
motion must also be two-dimensional. 

There is, however, an a priori dithculty with any theory which supposes 
that the perturbation remains small. Since the circulation round any 
material circular circuit concentric with the axis must remain constant, the 
radius of such a circuit must always be nearly equal to its initial radius. 
At first sight this restriction on the total strain of the fluid seems unlikely 
to be satisfied since it is inconsistent with any prolonged general streaming 
(however small) past the body. Nevertheless, Stewartson’s (1952, 1953) 
solutions do show that the perturbation can remain small, and these solutions 
must theiefore contain an explanation of the mechanism by which the 
circulation is maintained at a constant level, though Stewartson does not 
point this out. Moreover, his solution is very complicated, largely due to 
, and this rather 


Y 


the formation of the singular surface at the cylinder ( 
obscures the mechanism of the flow. In this note, therefore, a much simpler 
solution which does not have a surface corresponding to C is obtained, the 
primary aim being to illustrate the mechanism by which the circulation 
remains constant even when the perturbation remains small. 

The flow considered is that due to the slow uniform motion, started 
impulsively from relative rest, of a paraboloid of revolution along the axis 
of a rotating liquid. It is found that the radius of any material circuit 
concentric with the axis of rotation executes small oscillations which are 
180° out of phase with the corresponding oscillations in the azimuthal 
velocity component. As time progresses, the amplitude of these oscillations 
decays to zero everywhere except on the paraboloid. ‘The ultimate flow 
is then steady and two-dimensional everywhere except on the paraboloid 
and on the axis of rotation. On the paraboloid the velocity oscillates 
finitely and on the axis the velocity component parallel to the axis oscillates 
finitely and the other components are zero. ‘The swirling motion about the 
axis found in the case of the sphere and ellipsoid is absent here. 

It may be pointed out that the results obtained here can be deduced 
from Stewartson’s (1953) solution for an ellipsoid by carrying out the 
usual limiting process. But the procedure adopted here is found to be 
simpler. 

2. SOLUTION TO THE PROBLEM 

We choose cylindrical polar coordinates, Oz along the axis of rotation 
and (r, @) polar coordinates in a plane normal to Oz. Let the unperturbed 
motion of the fluid consist of a uniform angular velocity Q about the z-axis. 
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A paraboloid of revolution (whose axis of symmetry coincides with the 
z-axis) impulsively starts to move along the axis at ¢ = 0 with uniform 
velocity V. 
have in effect superposed a uniform velocity — V on the system and brought 
Let the components of the fluid velocity along the 


If we choose the origin of coordinates to be in the body, we 


the body to rest. 
directions of increasing 7, 0, 2, be u, Qr+v, w, respectively, where u, v, w 
are small. ‘Then the linearized equations of motion are 








u Or | 
ae —20v =- =~, 
ot or 
ov 5 5) 
ot T 20du = 0, f (2.1) 
ow oP 
Of -~ fakes m 
lo u 
-—(ru)+— = 0 2.2 
7 > ) Os s ( ) 
where p a 1 Q2;2 
; Rar. 


‘The boundary conditions are that 

u—> 0, vy —>0, w—>—V, as 2-—> 00 for fixedr,t, (2.3) 
and, on the body, the component of the fluid velocity normal to the body 
iS Zero. 

As Morgan (1953) pointed out, the initial disturbance travels with 
infinite velocity and the initial motion relative to the rotating system must 
be the irrotational motion with the given boundary conditions. ‘Taking 
the velocity potential of this irrotational flow to be 

f(r, 3) -Vz+x(r, 2), 


we have at ¢ ) 


Now to take the Laplace transforms of u, v, w, and P, we put 


u | u(r, 2,t) dt, etc 
Chen (2.1) and (2.2) become 
su—2Qv = — oA : 
or 
yi 2Ou Q, (2 4) 
oN 
STH ] = . 
lo Ou fe 
-~—(ru)+ — = 0, (2.9) 
r or Oz 
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where N = P-x; 
and the boundary conditions (2.3) become 
— - = I 
a 2 ee = as Z—> 00. 





n= om 402 7 ) 
I LON 
w=-—---e 
$ § OZ 





12 (aM) , eesanes 
= r—)4 ae 7a = U, 


ror 


r sz 
/ 
and the boundary condition (2.6) becomes 
ol if 
ON 50 asz>o. 
Os 


If the section of the paraboloid in the (z,r)-plane is 2 
condition on the body is 


2arut+w = 0, or 2Zaru+w = 0, 
or, in view of (2.7), 
sy ON. ON 
Zar =——s = - J 
S°-- #82 OF O02 


So equation (2.8) is to be solved with boundary conditions (2.9) anc 

Now we can easily formulate the problem in a coordinate sy 
which (2.8) can be solved simply and in which the body is a co 
surface, by taking a suitable transformation of independent v 


We introduce new coordinates (€, 7) defined by 


where 





On the paraboloid we have 
& = & = (K/4a)"2, 


With the above transformation equation (2.8) becomes 





o2N LoN , 2A idN 0 
o&* & ( e ( 7” ? 1) on : 
and the boundary condition (2.10) becomes 
oN — —_ 
—~- =—2KVE on € = Se 
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(2.6) 
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(2.10) 


1 (2.10). 
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iriables. 


(2.11) 
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The appropriate solution of (2.12) is 
N = (4+ Blog é)(C + Dlog yn), 


and using (2.13) we get 


N = —2KVé£, log é. (2.14) 
Now 
oN ! oN oN V &2 
ez = 2K(€?+ 7) (¢ Oe 7] *) E24 2° 
Therefore dN/dz > 0 as €-—> w%, in agreement with (2.9). ‘Thus (2.14) 


is the appropriate solution. 

The results for u, 7, w follow immediately from (2.7). Finally, inverting 
these Laplace transforms, we find that the velocity components at any 
point of the fluid are given by 


ae L Age)c2 + 402\ esl 
Ae | (1 (1 + 4az)s* + 4Q )s ds: 


Sarai. w s 
QV prt (1 + 4az)s? + 402?\ e* : 
eee | ep aS (2.15) 
4arzmi , OF w* ce 
V prim 24402 
w=—) —. ——. e* ds 
271 | w~ 


where y > 0 and 
w? = {(1+4az)? + 16a2r2}"2(s2 + 4002/2) 2(52 + 4Q2/2)H2 


» 14+4az+8a?r? — 8ar(a*r? + az)? 

‘ (1+ 4az)? + 16a?r? - (2.16) 
>»  1+4az+8a?r? + 8ar(a?r? + az)" 
ere (1+ 4az)* + 16a?r? 








These results represent a complete formal solution to the problem. For 
the present purpose of ascertaining the general features of the flow, however, 
it is only necessary to consider certain special cases of the formulae (2.15). 


3. GENERAL FEATURES OF THE FLOW 
On the surface of the paraboloid the integrals (2.15) simplify considerably 
and it is possible to evaluate them in the forms 








2Var 2Qt 
Oe 14 4ar © (1 +4a?r?)!/2? 
2Va 2Qt 4 
; (1 + 4a? ara S19 + 4a2r?)l 2? (3.1) 
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9 ¢ S 9.9 9° 
1+ 4a?r? (1 + 4a?r?)t2 
‘Thus we find that the motion never becomes steady on the paraboloid. 
More important, perhaps, these results show very simply the way in which 
the circulation round circular material circuits, concentric with the axis 


of rotation and lying on the paraboloid, remains constant. Since the radial 
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velocity u oscillates sinusoidally in time, the radius of such a circuit must 
oscillate in a similar way, and the primary rotation then makes an oscillatory 
contribution to the circulation round the circuit. ‘his must, in turn, be 
counterbalanced by an oscillatory contribution from the azimuthal pertur- 
bation velocity v, in accordance with (3.1). As far as the validity of the 
linearized analysis is concerned, the essential feature of this mechanism 
is that no fluid particle is displaced appreciably in a radial direction from 
its initial position. 

Similarly, on the axis of rotation we find 

; 4Vaz 2Qt 
Uu 0, o= 0, w= - oe cos (T + 4az)ie* 
Here again the oscillatory axial velocity implies that small material circuits 
surrounding the axis of rotation are never swept on to the surface of the 
paraboloid, thereby increasing their perimeter by a large factor; an essential 
result if the azimuthal perturbation velocity is to remain small. 

From these special cases it is reasonable to infer that the same oscillatory 
mechanism is responsible for maintaining the radial positions of all fluid 
particles, and this may be verified directly when the motion is approaching 
its ultimate form. ‘Thus, for large values of ¢t, the integrals in (2.15) may 
be evaluated by inserting cuts in the s-plane from s = +27Q/, and 
s = +27Q/, along lines on which the imaginary part of s is constant and the 
real part decreases. ‘he path of integration may now be replaced by a 
path round the infinite semicircle #{s; <0 and round the four cuts. For 
example, the contribution from the branch point s = 27Q/, to the integral 
in the first of the equations (2.15) is found to be 
] - + 4az)l? (3) 
© 1(P=ByR4OLyY? A 


for large ¢t. In this way we find te 


A 1—(1+4az)l? 
“nm — aml - (1+ 4a2)l; - sin(2Q/, t— 4) + 


(3.2) 














16(ar)?(a?r? + 1,(Q22,)" 
1—-(1+4az)i2 . 1 
1 Ql, ye sin(2Q/, t tia lo 
V 1—(1 + 4az)/? 





U~ ya cos(2Q/, t— 4) + 
16(ar)?(a?r? + az)"4 | I (a3h,)"? ennt~T) 
_1-( +4aale cos(2Q/, t + 477) 
13(QQ/,)"? —— a a 
V 1—/? 
L_ sin(2Q/, t — 42) + 
wn 75-5 5 sin(2Q m 
4(ar)'?(a*r? + az)" | ars pe 20) 
1--72 . l 
i(OL) sin(2Q/, t T 7) | (at) 
Thus the only significant difference here is that the amplitude of the 
y sj Pp 
oscillations decreases to zero, so that the ultimate motion is in general 
steady and two-dimensional and the axial velocity of the fluid is ultimately 


the same as that of the paraboloid. 
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In view of the ultimate singularity in the velocity gradients on the axis 
and body it seems that the detailed form (but probably not the general 
nature) of the solution is of doubtful validity in this neighbourhood. 

In conclusion, I wish to thank Professor B. R. Seth for his kind guidance 


throughout the preparation of this paper. 
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On the series that represent tides and surges 
In an estuary 
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SUMMARY 

This paper is concerned with a progressive wave of general 
form in an infinitely long estuary of uniform cross-section when 
there is a permanent current independent of the wave. ‘The only 
approximation is the neglect of friction. Explicit formulae in 
the form of infinite series are found for the surface-elevation and 
for the current. 

In the special case in which there is no permanent current and 
when the oscillation at the mouth of the estuary reduces to a single 
harmonic constituent, the first five harmonic shallow water 
constituents at any place up the estuary are evaluated. 


1. INTRODUCTION 


In the making of tide-tables for a port, the basis of the method usually 
employed is that of a series of harmonic constituents. When the port is 
on an estuary, a number of these constituents are of the type usually called 
‘shallow water’ constituents. For many years Doodson (1957) has used 
his method of ‘harmonic shallow water corrections’ in predicting the 
times and heights of high and low water. ‘This method implies the existence 
of a much larger number of shallow water constituents than the number 
of harmonic units on any predicting machine. 

It is the object of the present paper to examine the series in the case 
of an infinitely long estuary of uniform cross-section when there is a 
permanent current down the estuary and no friction. ‘This is an ideal 
problem in the three respects of infinity, uniformity and absence of friction. 
But, because precise results are obtained, the investigation appears to be 
of value in spite of the limitations indicated. 

As the cross-section of the estuary is uniform, the permanent current 
is also uniform. Owing to the absence of reflection, the motion takes the 
form of a progressive wave, and no serious limitation is placed on the shape 
of this wave at the mouth of the estuary. 

The simple basin specified is that which was taken by Saint-Venant 
(1871), McCowan (1892) and Fjeldstad (1941) when they evolved their 
accurate finite solutions of the fundamental equations. Saint-Venant did 
not include a permanent current, but McCowan did. ‘Their solutions show 
that the heights of high and low water at any place up the estuary are the 
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same as at the mouth, and this gives a complete answer to the most important 
question regarding a combination of tide and surge. ‘They also show that 
a bore will always form, the position of its formation depending on the size 
and shape of the oscillation at the mouth. 

But the solutions of Saint-Venant, McCowan and Fjeldstad do not 
express the surface-elevation as an explicit function of position and time, 
and it is just such a function which is required in tidal prediction and which 
is provided by series. 

I first obtained the general formulae of the present paper by substituting 
series with unknown coefficients into the fundamental differential equations, 
But Dr M. S. Longuet-Higgins then pointed out that the results could be 
obtained much more shortly from the implicit solutions by means of 
Lagrange’s theorem on implicit functions. I now follow this method. 

‘The early terms of the series have an interest in themselves, and they 
take simple forms when there is no permanent current. Taking the special 
case in which the oscillation at the mouth reduces to a single harmonic 
constituent, the first five harmonic shallow water constituents are evaluated 
for any place up the estuary. The first three of these were given by Airy 
in 1842. 


NOTATION AND GENERAL EQUATIONS 
Denote by : 


the acceleration of gravity ; 


4S 


h the undisturbed depth of water; 

v distance up the estuary from the mouth; 

t the time; 

the elevation of the water-surface above its undisturbed level ; 

uw the current up the estuary ; 

uy the permanent current down the estuary, assumed constant and 
uniform ; 


F(t) the ratio C/A at the mouth of the estuary; 


ind write Te 
( (gn), ] v/(C— Ug), n= t—@, 


"he equation of continuity 1s 


\ OG 1 
\/2 + O)U; = U (1) 
x t 
ind the equation 1 motion is 
Uu u ¢ 
) 
7 9 (). 2) 
XY x 
while a stationary solution of these equations 1s given by 
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Following Lamb (1932, §187), multiply (1) by [g/(h+{)]'* and 
subtract from (2); the resulting equation may be written 











ou g \l20¢ 

— - y yee ( ) + 
ot (; [e :) ot , (4) 
ou g \'20¢ 2 
Sarg ac : —— = ( iS 
ax (; + :) a= (>) 


Both (4) and (5) may be integrated to give 
U+U, = 2[g(h+C)]!? —2c, (6) 


where the constant of integration has been chosen so as to satisfy (3). 
Equation (6) may be written in the forms 





U+Uy C\ 12 ss 
eh | eee ae / 
c ( i) , 7) 
C utu 1 /ut+u,\? 
$ _ Pas 0)" (8) 
h c 4\ ¢ 
Substitution from (8) into (2) gives 
ou a e.. ' 
7 +(c¢+ uy + 2u) = ze 


and into (1) gives 


u+Uu,\ (ou ou 
| 4: Seger eat c+Aiu + 34 —+=(Q, 
( ac | ( 249 2 Vax | 


ou Ou 
y= als os U Yeas Q, 9 
ot OX ( 
where 
U = c—uy+ 3(ut+ up). (10) 


Substitution from (4), (5) for du/dt, du/dx into (9) gives 


at a 
—~ +U—=( 
> l x ), (11) 
and from (7) into (10) gives 
l C\H2 tty 
bash 2 ss 2 


The velocity U is that of propagation of an element of the wave up the 
estuary. 
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It is easily verified that solutions of the differential equations (11) and (9) 

are given respectively by 
C/h F(t—x/U), (13) 
(u+uy)/c = f(t—x/U), (14) 
where F(), f(.) are any differentiable functions. Since x= 0 at the 
mouth of the estuary, the function F is that which has been defined 


above. ‘The function F(t) will be regarded as prescribed, and from its 
physical meaning it follows that F(t) > —1. From (7) it follows that 
tf= (14+ F)*-1, (15) 


the arguments of F, f being the same. 

Equation (13) is equivalent to McCowan’s solution (1892), and when 
Uy = 0 it is equivalent to Saint-Venant’s solution (1871). 

It is seen from (12) and (13) that, for the same time, the ettect of the 
permanent current is to make a specified surface-elevation occur at a 
shorter distance up the estuary than when there is no such current. 


3. "THE SERIES 

‘The equations (13) and (14) give ¢ and w as implicit functions of ¢ and x, 
but it is the object of this paper to obtain ¢ and w as explicit functions. 
Such explicit functions, in the form of infinite series, may be derived 
from (13) and (14) by means of Lagrange’s theorem on the expansion of 
a function which is defined implicitly (Edwards 1896, p. 451). 

Lagrange’s theorem may be stated as follows. Let 4(y), ¢() be functions 
which may be expanded in powers of y — », and let 


y=nt0 (Vy) ; 
gn 
then b( y) = (7) 4 > =F [{h(7) |" PO(n) | 1) (16) 
n l 2 


where the bracketed indices denote derivatives with respect to 7. 
Now, in the present application of the theorem, the parameters 7 and 6 
are given the meanings indicated at the beginning of § 2, so that 


x C—o 
t —_ —_ = n + a(1 a ‘) = Y, 
l ~ 





U 
on taking dy) = 1- c = (17) 
thus (13) and (14) may be written as 
tih = Fly), (18) 
(u+uy)/c = f(y). (19) 


From (12), (17) and (18) it follows that 








is, (20) 


and from (15) that 

















(i+ “FC yy? 
it will be supposed that Fy) and f(y) satisfy the condition imp osed on dy) 
in Lagrange’s theorem. From (18) and (20), the expansion (16), with 
F in place of ys, becomes 





C J a , 
h F(n) | > f Z,(), (22) 
t n=1 
where 
L | [1+ F(y)}?-1 | 2 i 
Z,, (1 =—| <= ——— a > Fa (n (23 
wn) ail iii + F’(n)]!? — fu/e— §) | or] 


for n = 1, 2, 3, ..., and from (15), (19), (20) and (21), the expansion (16), 


7 = 9 


ith fi lace f &. beco 2g 
With 7 in Pp ace OF yw, Decomes 











— 2[1+ F(y)?-2+ > OU,(y), (24) 
where 
\ [14 < yr" — gUg/C— 3) (1 yn 
‘lhe series es ) and (24), with their coefficients given by (23) and (25), 


constitute the bs aria of the fundamental equations in terms of 
F(), which itself gives the prescribed elevation at the mouth of the 
estuary. 

The condition imposed on ¢(y) in Lagrange’s theorem requires that 
¢(¥) shall not be infinite, and, by (17), this requires that U shall not vanish. 
Assuming that U, the velocity of propagation up the estuary, is once 
positive, it follows that, for the validity of the expansions (22) and (24), 
it must always be positive. From (12) and (18) this requires that 

3[1+ F(y)]!?—uy/e—2 > 0. (26) 

When the expressions (23) and (25) for the coefficients Z,,(7) and U,,(») 

are expanded in ascending powers of F(y), the first term in each is 


Rr 3F(n) )n mae (n—1) 
ail 130 aot f (| . 


| | 3 \" 





and this is equal to 





(n+ 1! 20 - Uyic){ LAG)? 


Examples of this will occur in § 4. 


| The expansion of (1+ )'? in ascending powers of F requires that 
\F,\ <1. On combining this with the condition (26) it follows that, for 
the validity of the expansions, 
| Ug 2 
Af. ; y @ | 
Ab +2) 1+ F(n) <2. (27) 


\t the mouth of the estuary, the relationships (27) give the lower and 
upper bounds of the ratio of the total disturbed depth of water to the 
undisturbed depth. 
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4, EARLY TERMS OF THE SERIES WHEN THERE IS NO PERMANENT CURRENT 
When uw, = 0, then 





0 = alc, n = t-2x/€¢, 

and the early terms of the series (22) and (24) may be written as 
F 2 rd v4 
¢ ee ee: ee eee 
h = I + 74+ Bata 437 aly (28) 
u x x os a Q2 
- = Uy+ -U +—U,+—,U,+4 U 9) 
c oo ek ee oe ee” Oe 


/ 


respectively, where U, = 2(1+ F)!?—2, and the arguments of P; a, W. 


are t—x/c. 
On using (23) and (25) and expanding Z,, U,, in ascending powers of f 


as far as F°, the following results are obtained: 


3 oy 75 O77 (ay 
by = E F2~ * ee. #9 Fie 1077 F* : 








] 
4 8 64 640 | 
ee 1341 Je | 
Z, = E F3— 5 Ft+ ie ~~ | 
n(n oD. 
i oe E 3 |” 
640 ‘ 
nina P+ iP a Ft La oe (30) 
~ E — oem a re aa el, 
U, = E F3— : Fly = ey’, | 
: 207 2) | 
all | sa? 3a0 | > | 





27 a) | 

U, = ——frs| . | 

: 640 J 

These results may also be obtained from the fundamental equations by 
the method of successive approximations. 


5. HARMONIC CONSTITUENTS 
When there is no surge, the surface-elevation at the mouth of the 
estuary is given by a series of the form 


F(t) = > A, cos(o,, t+«,,), (31) 
a= | 


7 


where A,,, o,, %, are constants. When (31) is substituted into the formulae 
of §3, the surface-elevation at any place up the estuary is given by a triple 
series of harmonic constituents. 

When the series (31) reduces to one term, so that F(t) = A cos(ct +), 
then F(t—x/c) = Acos@, where 6 = o(t—x/c)+%. Substitution into the 
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results = §4 then gives 





Z, 21 ’ 

thas ~ | Atsin204 33 4° (sin 6 + sin 30) — ! 

75 ital ie ean oe | 

— 738 A*(2 sin 20 + sin 40) + 7048 A*(2 sin 9+ 3 sin 34+ sin 58), | 

La a 63 

ed 39 AX(cos 6+ 3 cos 34) + 35 A'(cos 20 + cos 48) — | 
$ (32) 

>S17 A5(2.cos8-+9.cos 38 50), | 

— 3048 (2 cos cos 34 +5 cos 58), | 

Bs 189 

ile s A*(sin 24 + 2 sin 44) — jo24 A°(2 sin 8 + 27 sin 36+ 25 sin 58), | 

oe val A*(2 cos 8+ 81 cos 36 + 125 cos 5@) ! 

ms ~ 2048 d ~ . 4 oe a oe . | 


The terms of the above formulae in A* and A® are equivalent to formulae 
given by Airy in (1842). 

Rearrangement of the series (28) so as to group together terms with the 
same function of @ gives 


C { 9 ax\2 1 ox\2) 1) 
= a — 422149~3(—) s | scosd+ 
Z=A\l- yA (=) | 1+ 354 14 3 (=) } |p eose 


toa 3 a safe LLs 
+ 35 AB 121 + 35 AP 1077-378 (=) |}sine+ 


r 
63 ox\2 
Seay 4) pes Ii = 
34 a cos 24 
3 ox ( 1 25 ox\2 ]} 
— — A*— 2] — —3( — 26 — 
44 ih ake 4 (=) |} sinzo 


4 J 

AS S) 41+ 4 149-27(2)' || cos30+ 
ae. 359 ox\27) . 

AB 374 aell = S61) tae 
Cc L 32 p c f 


ox\: 3 x” (25 x\2)  . 
+54(=) eth. = re -3(2) + sin 40 — 
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Some properties of boundary layer flow during the 
transition from laminar to turbulent motion 
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SUMMARY 

‘Lransition in the boundary layer on a flat plate is examined 
from the point of view of intermittent production of turbulent 
spots. On the hypothesis of localized laminar breakdown, for 
which there is some experimental evidence, Emmons’ probability 
calculations can be extended to explain the observed statistical 
similarity of transition regions. Application of these ideas allows 
detailed calculations of the boundary layer parameters including 
mean velocity profiles and skin friction during transition. ‘The 
mean velocity profiles belong to a universal one-parameter family 
with the intermittency factor as the parameter. From an examina- 
tion of experimental data the probable existence of a relation 
between the transition Reynolds number and the rate of production 
of the turbulent spots is deduced. A simple new technique for 
the measurement of the intermittency factor by a Pitot tube is 
reported. 


1. INTRODUCTION 


The stability of laminar shear flows and the transition to turbulent 


motion has received and continues to attract much attention because of 


its fundamental importance to the study of fluid motions. A great deal 
of theoretical and experimental work has been done on the instability of 
laminar boundary layers and in determining criteria for transition. It is 
known that in general transition can be influenced by the following factors : 
surface roughness; free stream turbulence; surface curvature; pressure 
gradient; surface temperature; Reynolds number; Mach number; 
acoustic radiation; injection or suction of fluid at the wall. Because of 
the complex manner in which the various factors affect transition, no 
satisfactory theory for the transition process has been possible so far, and 
the origin of turbulence is still largely an unsolved problem in fluid 
mechanics. 

Recently, recognition of the intermittent character of laminar break- 
down in a boundary layer has directed attention towards the fundamental 


processes involved. ‘lhe observations and theoretical work of Emmons 


(1951) and the experiments of Mitchener (1954), Schubauer & Klebanoft 
(1955), ‘Tani & Hama (1953), Hama et al. (1957) and others have shown that 
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the transition phenomenon in a boundary layer is characterized by the inter- 
mittent appearance of turbulent spots which move downstream with the fluid. 
The mechanics of spot creation and growth is as yet not completely clear. 
Experimental observations with flow visualization techniques suggest that 
the amplification of ‘Tollmein—Schlichting waves becomes associated at 
some stage with the concentration of vorticity along discrete lines, which 
subsequently distort into vortex loops in the boundary layer. ‘The vortex 
loops themselves go through a process of distortion and extension finally 
resulting in the creation of ‘spots’ of turbulence. Once created these 
spots are swept along with the mean flow, growing laterally as well as 
axially with laminar flow in their trail. he spots originate in a more or 
less random fashion and increasingly overlap as they enlarge during their 
transit downstream, finally covering the entire plate and resulting in fully 
turbulent motion. Passage of the spots over points on the surface results 
in alternations of laminar and turbulent flow. 

These alternations can be quantitatively described by an intermittency 
factor y which represents the fraction of time any point spends in turbulent 
flow. Flow at zero pressure gradient over a flat plate is the only case so far 
studied in any detail. For this case, when transition occurs naturally or 
due to a disturbing agency which causes the spots to appear at some distance 
downstream, the spots grow laterally in a nearly linear manner sweeping 
‘turbulence wedges’ on the plate. ‘The angle of these wedges is approxi- 
mately the same as that for transverse contamination found by Charters 
(1943). During the initial period the spot growth is non-linear, and 
envelopes of spot growth exhibit a characteristic curved shape. 

The existence of turbulent spots in boundary layer flow has been 
confirmed experimentally by Mitchener (1954), Schubauer & Klebanoft 
(1955) and in this laboratory, and it appears to be quite certain that their 
occurrence plays a fundamental role in the mechanics of boundary layer 
transition and probably also in the breakdown of laminar motion in general. 
In this paper the transition region of a boundary layer on a flat plate is 
examined from a viewpoint which keeps the intermittency of flow as the 
central feature of the process. It is then sought to utilize the available 
experimental and analytical information to explain details of the flow in 
the boundary layer during transition. ‘The main sources of experimental 
information on the mechanics of transition utilized for analysis here are the 
experiments by Schubauer & Klebanoff (1955) at the Bureau of Standards, 
Washington, D.C., and those conducted at the Indian Institute of Science, 
Bangalore. Both sets were carried out in low-turbulence wind tunnels 
under carefully controlled environmental conditions. 

In the literature it is indeed rare to find transition measurements which 
provide sufhciently complete information. ‘The most serious omission in 
the older studies is usually the intermittency factor whose significance was 
probably not recognized earlier. It is hoped that such measurements will 
become more common in the future, and will provide more extensive 


data on transition flows. 
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2. DISTRIBUTION OF INTERMITTENCY DURING TRANSITION 

Experimentally, distributions of the intermittency factor during a 
transition zone can be measured by placing a hot-wire close to the surface 
and recording the amplified output on an oscillograph. ‘Typical records 
obtained in this manner are shown in figure | (plate 1). The passage of each 
spot over the hot-wire is characterized on the oscillographs by random 
high-frequency fluctuations in contrast to the relatively regular, smooth 
laminar regions. ‘The factor y is easily calculated from such records by 
counting the intervals of time for which the traces show turbulent flow. 

For a quantitative description of the transition process, Emmons (1951) 
has formulated a theory based on probability considerations. This theory 
assumes the existence of a turbulence spot production function g(x, y, t), 
and Emmons shows how this can be related to the probability of the flow 
being turbulent at some point, i.e. to the intermittency factor y. It has 
been shown by Narasimha (1957) how the original assumptions about the 
production of turbulence sources have to be modified in the light of new 
observations. ‘The main feature of the modification stems from observations 
which suggest that laminar breakdowns in a two-dimensional flat plate 
boundary layer are very nearly point-like, and that the spots originate in 
only a restricted region. On physical grounds which have some support 
from stability theory, one may argue that the amplified laminar oscillations 
can be expected to lead to production of the turbulent spots at a fairly 
definite distance from the leading edge of the plate. Since the actual 
occurrence of spots is random, one can reason that there must be a region 
on the plate downstream of the point of instability where the origin of the 
spots is most probable. At points further downstream, the tendency for 
breakdown rapidly decreases on account of the state of flow behind a spot 
as well as the fact that most of the breakdowns have already occurred. 
Thus the source-rate density function describing spot production should 
have a maximum at some location. It turns out that the experimentally 
observed beginning of transition marks this location. 

If one assumes a Gaussian error curve with its maximum at x; to describe 
the turbulence source-rate function, then Emmons’ theory can be used to 
calculate the resulting y distributions. Here x,is the location where transition 
begins. It is interesting to note that the Gaussian curve which gives the 
best agreement with experimental data has a standard deviation approaching 
zero. ‘This implies the form of the source-rate function to be very nearly 
the Dirac deltafunction. Figures 2 (a) and 2 (6) show representative cases of 
transition illustrating calculations based on the above arguments. An 
interesting and important feature of this hypothesis of localized laminar 
breakdown is the explanation it suggests for the observed statistical similarity 
in the transition distributions. Application of this hypothesis together with 
Emmons’ theory gives 


y= 1-e%*. (1) 


Here € = (x—x,)/A, a normalized stream-wise coordinate in the transition 
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zone, A being a measure of the extent of the transition region characterized 
1 by 
: A = (X}, 20.75 — (fy = 0-259 
S and A is a constant equal to 0-412. 
L Equation (1) shows that whatever the cause of transition, all transition 
regions define a single universal intermittency distribution. Figure 3 
L shows the available data on intermittency during transition, and demonstrates 


the agreement between equation (1) and experiment. Further experimental 
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Figure 2. Calculated y-distributions for normally distributed source-rate density. 


evidence supporting the view that the origin of most of the turbulent spots 
is very close to the beginning of transition is given by observations of 
boundary layer thickness during transition. Figure 4 shows two repre- 
sentative cases. They reveal that the virtual origin of the resulting turbulent 
boundary layer is very nearly at x, The effect of pressure gradients on the 
y(é) distributions has not yet been investigated in detail, but preliminary 
measurements mdicate that the distribution is modified mainly in the 
neighbourhood of y = 0. Pressure gradients would thus appear to chiefly 
influence the manner of spot growth in the initial period. 
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Another point of interest is the distribution of y across the boundary 
layer, viz. normal to the plate surface. Note that y(y) is a function of the 
shape of the turbulent spots. ‘The very detailed experiments of Schubauer 
& Klebanoff (1955) and measurements in this laboratory show that the spots 
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Figure 3. Universal distribution of y vs € with transition due to different agents 
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Figure 4. Boundary layer thickness in transition region (arrows mark limit of 
logarithmic turbulent boundary layer). 


have a nearly constant cross-sectional area close to the surface, but they 
taper towards the outer edge. ‘The measured y(¥) distributions + 
transition are similar to those observed by Corrsin & Kistler (1954) i 

fully developed turbulent boundary layers, i.e. y varies from a constant 
maximum value close to the wall to zero towards the edge. While this 
y(y) variation is probably of importance to the detailed structure of the 
turbulent motion associated with the spots for all practical purposes, the 
y value near the wall is the characteristic property of importance for the 
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transition region. For instance, it is seen in §3 that the y(y) variation has 
: only a secondary influence in determining the mean velocity profiles in 
transition. 

‘lo sum up, it would appear that the almost point-wise breakdown of 
laminar flow in a boundary layer while occurring randomly in time takes 
place very nearly on a discontinuous line across the flow, whatever the cause 
of breakdown may be. 


) 


3. DETAILS OF TRANSITION FLOW 


(a) Location and extent of the transition zone 





It has been observed in §2 that the distribution of intermittency is 
universal on the y(€) plot, ‘his implies that the transition zone has always 
the same extent in terms of the €-coordinate. ‘This information tells us 
nothing about the physical location of the beginning of the transition region 
So far it is not possible to predict theoretically, even in the simplest case, 
the Reynolds number of transition R, Since the mechanics of creation 
of the turbulent spots has yet to be formulated, attempts such as those by 
Liepmann (1945) and Lees (1952) to estimate R, from the development of 
amplified laminar oscillations lead to difficulties and have not yet been 
too successful. Experimental investigations have established useful 
correlations for the effect of particular types of disturbance agencies on R, 
(see Dryden (1953), Gazley (1953) and Lin (1955) for recent reviews of 
the available information). ‘There are, however, unexplained discrepancies, 
and in any case these correlations are not general enough to predict R, 
when there is a combination of transition agencies. 

A recent study of the experimental data on boundary layer transition 
available in the literature shows the possible existence of an interesting 
relation between the transition Reynolds number R, and R,, the Reynolds 
number based on the physical extent of the transition zone. Figure 5 
shows a plot of R, vs R, for the available experimental data. Since in 
many cases there were unspecified differences in experimental conditions, 





the data show considerable scatter; but in spite of this one can discern 
the dependence of R, on R,. ‘The extent of the transition zone appears to 
decrease progressively with decreasing R,, to some limiting value. Actually 
one would expect a family of R; vs R, curves with each curve depending 
on the specific agency or agencies causing transition and each having perhaps 
different limiting values. Figure 5, however, shows that this dependency 
on the type of disturbance is not too marked, and that for rough estimates 
a single mean curve may be used. ‘The data are represented on the average 
by the equation 


R, = aR? (2) 


with a = 5-0 and B = O'S. 





It should: be noted that R, is related to the turbulence source-rate 


production function of Emmons, and the possible existence of an R, ws R, 
relation is probably of some significance to the transition problem. Using 
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equation (2), it is possible to calculate %, the average number of turbulence 
sources produced per second, as a function of the transition Reynolds 
number. The results of this calculation are shown as an #-scale in figure 5. 
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Figure 5. Correlation between R, and R, (Flagged points are 
hot-wire measurements.) 


Subsonic measurements Supersonic measurements 
@ Indian Inst. Sci. 7 Coles (1954), flat plate, 
Schubauer & Skramstad | M 1 -97-4-54, 
(1947), Schubauer & | (| Higgins & Pappas (1951), flat 
Klebanoff (1955). > flat plate plate, M = 2-4. 
4. Dhawan (1953). / Chapman & Kester (1953), 
’ Hunsaker (1939). cone-cylinder, M 2-9. 
Burgers (1925). J (Winter et al. (1954), cone, 
@ Wright & Bailey (1939), ME =-2:4553°25. 
flat plate with pressure gra- 
dient. 
‘Y Silverstein & Becker (1939), 
aerofoil. 


(b) Mean velocity profiles in the boundary layer during transition 
Following the discussions of §1 and § 2, we note that each point on the 
flat plate during transition is covered with a laminar boundary layer except 
during those intervals of time when turbulent spots pass over it. ‘T'wo 
questions of immediate interest are: (i) what velocity profile is associated 
with a turbulent spot during its motion; and (ii) whether the mean velocity 
profiles during the intermittent transition process exhibit any similarity. 
The clue to these and other related questions is contained in the 
observations which show that the turbulent spots are randomly created 
at a fairly definite x-wise location, and are on the average all similar in shape, 
their propagation cones being very nearly true cones with straight generators. 
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The fully turbulent boundary layer on the plate which results after transition 
exhibits the characteristics of a two-dimensional turbulent layer, because 
it is a random mix-up of the turbulent spots originating from a narrow 
localized region. Within the transition zone one can similarly associate, 
on the average, a two-dimensional boundary layer with each spot, since 
the transition boundary layer is also a random mix-up of the spots but for 
smaller intervais of time. ‘The overall effect, therefore, is approximately 
equivalent to a flow in which, at each point, the boundary layer alternates 
between a two-dimensional laminar and a two-dimensional turbulent layer, 
each at its appropriate Reynolds number, y representing the proportion 
of time spent in turbulent motion. ‘The origin of the laminar layer is at 
the leading edge of the plate, while y departs from its zero value at some 
point further downstream. ‘The thickness of the transition boundary 
layer would be equal to the laminar or the turbulent layer starting from ~,, 
whichever is the greater. 
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Figure 6. Boundary layer velocity profiles in transition. 


Analysis of carefully obtained experimental data confirms the above 
observations. Figure 4 shows two representative cases of transition flow, 
and it is seen that the observed thickness of the transition boundary layer 
closely checks estimates based on the above reasoning. Figure 6 shows 
a case of transition analysed in detail. ‘The calculated laminar and turbulent 
profiles between which the flow alternates are shown for each station on 
the plate where a mean profile was measured with a total-head tube. Also 
shown on the figure are the calculated mean Pitot profiles obtained by 
mixing the laminar and turbulent flows. 

The agreement between the measurements and the calculations supports 
the basis of the arguments outlined above. ‘The calculated turbulent 
boundary layer profiles are based on a combination of the well-known 
‘law of the wall’ and the ‘velocity defect law’ as discussed by Clauser 
(1956). The constants entering into these laws were determined in each 
case from the measured fully developed turbulent profiles. Figures 4 and 6 
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show that the logarithmic turbulent profile is attained in the turbulent spots 
quite early in the transition process. ‘The behaviour of turbulent boundary 
layers at low Reynolds numbers being not very well understood, the deviation 
of the calculated profiles from the experimental data at low y values is not 
surprising. ‘lhe calculations can give only an approximate indication of 
the state of flow during the early history of a spot. 

The Pitot profiles u,(y) during transition were calculated from the 
equation 

up(y) = {1 —y)uit yup}, (3) 

y being measured independently by hot-wire techniques and u,, u, being 
the theoretical laminar and turbulent profiles respectively. Equation (3), 
incidentally, indicates that in regions of intermittent flow care must be 
exercised in interpreting the readings of averaging instruments such as 
Pitot tubes and hot-wire anemometers. ‘The reading of a Pitot is pro- 
portional to pressure (u?); hence the average velocity during transition 
flow indicated by a Pitot tube (i.e. equation (3)) is different from the true 
time mean #, which is given by 


c= —y)u;+ yup. (+) 





Equations (3) and (4) are compared in figure 7, which shows a representative 
transition velocity profile measured carefully by means of a small flat 
total-head tube. A complete y(y) distribution measured across the boundary | 
layer at the same location is also included. 
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Figure 7. Transition profile. 


Significant facts which emerge from the data and calculations of 

figures 6 and 7 are as follows: 

(a) The calculated Pitot profiles are in remarkable agreement with the 
reading of the total-head tube. This agreement proves the soundness 
of the physical picture on which the analysis is based. 

(b) Although y varies across the boundary layer, for purposes of profile 
calculation, the value measured close to the wall gives sufficiently 

{ for the whole profile. Thus the spot shape which is 


1 
e results 


accura 


associated with the y(v) variation appears to have onlv a secondary 


influence on the transition flow, y(x) being the significant property 
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(c) ‘Transition velocity profiles measured with a total-head tube must be 
corrected to give the mean profile. The errors are specially significant 
close to the wall, and skin friction determinations from uncorrected 
measured profiles are likely to have large errors. 

Before taking up the question of similarity of the mean velocity profiles 
during transition, it is useful first to consider the boundary layer parameters 
5*, 6, H which are known to play an important part in the quantitative 
description of boundary layers. During transition the displacement and 
momentum thickness can be defined as 


, | °1 (1 Uu ; (5 
oO, | nt Oe’ )) 
Ik Bee 
and “or ua / EN, 
0, = | U) dy, (60) 


where 6, is the thickness of the transition boundary layer, and ay) is given 
by (4). U’ is the uniform velocity outside the boundary layer 

Using (4) in (5) and (6), and recognizing that 6, = 6, if 6, > 4,7, or 
oO, 


5, if 5, > 5,, we obtain after some simplification 


6, = (1 yo, ty(d7), (3 a) 
and 
t= 1—y){(1—y)0,— yor} +ytver— (1 — yon} + 2v1 — y) FG), (6a) 
where df, 57 are displacement thicknesses and 6,, 6, are momentum 
thicknesses of the laminar and turbulent profiles respectively, and 


ihe Uu Uu 
F(8,) = 1 facbie ly. 
0 iL (a) (3 42 


Figure 8 shows two cases of transition flow for which equations (5a) 
and (6a) have been compared with experimental data obtained from 
corrected profiles and y measurements. One set of corrected profiles 
measured for R, = 4:3 x 10° is shown in figure 9, while the other set for 
R, = 2:3 x 108 is the same as that shown in figure 6. ‘The analysis is seen 
to give good agreement with the measured distributions of 5, and @,. The 
profile-shape parameter H is given by 


H, = 5, /6,, (7) 


and this can be put in the form 





A, r ) —, (7a) 
or 2y(1 
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where H, and H, are the form parameters for the laminar and turbulent 


boundary layers respectively 

ne variation of Fi is shown in figure 6 for the WO CdSCS studied. lt ] 
interesting to note that when plotted against the €-coordinate, the variation 
of H, in the transition region shows a near universality Actually, as th 
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figure shows, one can notice a small Reynolds number effect on H, mainly 
towards the latter half of transition. This is to be expected since the 
turbulent boundary layer associated with the spots depends on RK,, the 
Reynolds number of transition, and H for the turbulent layer is known to 
vary with Reynolds number. Figure 8 shows that the effect on H, is small 
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Figure 8. Boundary layer parameters during transition. 
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Figure 9. Corrected mean transition profiles (Indian Inst. Sci.). 


for a wide range of R, values. For approximate calculations, a simple 
expression can be derived for the distribution of the transition-shape 
parameter. This is of the form 


H, = Hp+(H,— Hye. (7b) 


‘This expression implies a linear relation between H, and y, and is plotted 
in figure 8 for comparison. 
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Having seen that a knowledge of R, and y(x) for any case of transition 
leads to a calculation of the mean velocity profiles as well as the parameters 
5*, 6, and H;, we now attempt to represent the profiles as far as possible 
ina unique manner. It is known from boundary layer theory that for steady 
flow past a flat plate with zero pressure gradient, both the laminar as well 
as the fully developed turbulent boundary layers exhibit similarity. ‘This 
fact is usually expressed by the profiles being independent of Reynolds 


number effects when expressed in the form 


u Yo , 
(7c), -() " 
u y \ ‘i 
and (3), =< (z)- (9) 


The function f() is defined by the Blasius solution, while g( ) can be 
tabulated in terms of the nearly universal turbulent boundary layer 
logarithmic laws, of which the !th power law is an approximate repre- 
sentation. Since the mean transition profiles are a resultant of the 














Figure 10. Similarity of transition profiles. 


laminar and turbulent profiles combined in proportion to y, it is logical 
to extend the similarity of equations (8) and (9) to transition flow by 
including the effect of y. We thus postulate that 


u - 
(7), = a7): a 


In principle the function A( ) is defined by a solution of the differential 
equations governing the mean transition flow. Unfortunately, these have 
yet to be formulated. We resort to an examination of experimental data 
for information on the form of A( ). 

In figure 10 two sets of mean transition profiles measured for widely 
differing R, have been plotted in the form (#/U), vs € for fixed values of 
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y/@. Since € is a function of y, a y-scale is also shown. Figure 10 shows 
that the transition profiles do indeed define, at least approximately, a 
function h(y/, y) which for all practical purposes is independent of Reynolds 
number. Figure 11 is a cross-plot from figure 10, and shows the universal 
laminar and turbulent transition profiles as well as the one-parameter 
family of profiles. Persh (1955) has noted, from an examination of experi- 
mental data, that transition profiles exhibit approximate parametric similarity 
when plotted as u/U vs H for fixed y/@. We recall from (7b) that H,(é) 
is nearly independent of Reynolds number effects. Since the functions 
y(€) and h(y/4,y) are universal, the reasons for the correlation observed 
by Persh become clear. The representation of transition profile data 
in the form of equation (10) is to be preferred as it shows the dependence 
on intermittency explicitly. 





a a 








Figure 11. Similarity of transition profiles; (u/U)), h(y/0,, y). 


(c) Skin friction 
Local skin friction 

If Cy, and Cry are the local skin friction coefficients for laminar and 
turbulent flow respectively, then the mean local skin friction during 
transition is given by 

Cy, = (L—-y)ey, + ep (11) 

This equation is plotted in figure 12 for the two cases of transition analysed 
earlier. ‘The transition c;-curves, which fair smoothly into those for laminar 
and turbulent flow, show a more rapid increase of wall friction than is usually 
assumed. At the end of transition the local skin friction is seen to be 
considerably higher than the values given by the turbulent skin friction 
laws based on the turbulent boundary layer starting at the leading edge 
of the plate. The usual presentation of skin friction data on a plot of the 
total friction coefficient C, vs R, obscures the fact that the origin of the 
turbulent boundary layer is approximately coincident with the beginning 
of transition and the onset of intermittency. ‘The usual skin friction curve 
for turbulent flow must be translated along R,, so that its virtual origin 
coincides with the foot of transition. 
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Independent evidence supporting this is provided by the transition 
local skin friction data of Coles (1954), which is shown in figure 13. Coles’ 
experiments measured cy directly, and were conducted at supersonic 
speeds so that the final turbulent boundary layer curves must include 
compressibility corrections. In figure 13 the Frankl-Voishel (1943) 
theoretical curve for c;,, has been used as this is known to give reasonably 


accurate predictions of supersonic skin friction. The characteristic shape 
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Figure 13. Local skin friction during transition in supersonic flow (from Coles 1954). 


of the local c; curves rising above the turbulent friction curves is seen again. 
These data also suggest that the mechanism of transition is not altered at 
supersonic speeds. The fact that after the breakdown of laminar flow 
the resulting turbulent boundary layer has an effective origin close to the 
beginning of transition resulting in high c, values is probably also the 
reason why measurements of transition recovery factors show values higher 
than expected (Shoulberg, Hill & Rivas 1954). 
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Another interesting feature of the transition local skin friction curves 
is brought out by figure 14, which shows that the difference {cy — Cy, =o} 
is independent of Reynolds effects during transition. Here Crly =o 18 the 
value of c,, at x, The points shown on figure 14 were calculated from the 
corrected experimental profiles used elsewhere in this paper. They 
approximately define a single curve with € as the variable. While it is 
dificult to give an a priori theoretical justification for this observed 
correlation, a reasonable explanation is as follows. 
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Figure 14. Similarity in variation of skin friction during transition. 


We have that in (11) c,, and c;,, may be represented as 
cr, = K,/Ry”, (12) 
Cr, = KA/RY, (13) 
where K,, A, are constants. Using the definition of € and equations (1), 
(12) and (13), we can write (11) as 
Cy, = FLR, R,, §}- (14) 
Now if R, and R, are related as suggested by the experimental data of 
figure 5, equation (14) becomes 
Ch = FUR, §}. (14a) 
The correlation observed in figure 14 shows this function to have the 
form 


FER, £} = O(R) + ¥ (14b) 
with @(R,) = cz, 9, the value of laminar skin friction at the beginning 
of transition. 
We thus get 
CA— fly =0 = ¥(E), (15) 


¢ 


function of the €-coordinate alone. We note that figure 14 implies, once 
again, that the main Reynolds number effect on transition is through the 


location of the region of laminar breakdowns. Subsequent events are 
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controlled by the intermittency and its distribution. An approximate 
expression for the function ‘YY is 
¥'(€) = 0-003(1 —e-**). (15a) 


This equation is plotted in figure 14. 


Total skin friction 

The behaviour of the local skin friction coefficient during transition 
discussed above can be used for engineering purposes to calculate the total 
skin friction of a plate on which transition occurs at some location. 
Integration of (11) for a plate of length / with transition starting at x = x, 
and ending at x = x, gives 


1 rl } 
CC. ==) ¢, dx 
; Ll! 0 
x, (2 A ie : , 
/() “0 
a (ee 1 , , 
~ T |p, det 1 | Cy dE + T | Cty 4%, for! > x. 


/ 0 J 


Here the upper limit of the second integral for / > x, has been taken to 
be € = 4-0 in accordance with figure 3. 

These expressions may be put in a more convenient form by utilizing 
the approximate relation (15a) and the R,vs R, relation of (2). The 
total skin friction coefficient C, is then given by 
: a+bR—C erf(dR’'/R’) 

ied R, 
and (16a) 
R'Cy,—k 


C, = ——>—_ for! > x, | 
: R, "4 


for! <x, 





where 
a= RiCp,\, =o} b = 0-003 + Cp, |. -0 


C= 0-00376%R/? = ()-0188R)”, d = 0:707/a = 0-141, 
b = C—a + 4aRi{Cp,, |x; —b} = 5{0-00376 + 4(Cp, |x; — 6) }Ry*— a, 


erfz = —| e * ds, 


Equation (16a) is plotted for a range of R, values in figure 15. For one 
case a comparison is shown with Emmons’ calculations, the usual 
Prandtl-Schlichting estimation, and Geber’s experimental data (see 
Schlichting 1955). ‘The present analysis shows a sharper rise of Cy, due to 
transition than is indicated by Emmons and shows a rounding-off of the 
transition curve at the beginning of transition unlike the Prandtl-Schlichting 
semi-empirical curve. 
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Figure 15. Mean friction coefficient of a flat plate with transition at R,. 


4+, ESTIMATION OF y BY A PITOT TUBE 
It has been shown in § 3 that fairly good agreement is obtained between 
the velocity profiles measured by a Pitot tube and those calculated from 
the formula (obtained immediately from (3)) 


(z) - lc (a). (ZF 
so - Of}, je! 


using measured values for y and the theoretically known values for (u/U), 
and (u/U),. ‘This immediately suggests the possibility of experimentally 
obtaining y without having to use the hot wire and associated equipment. 
From the measured velocity profile, or indeed from a single measurement 
of (u/U), in the transition region at a certain height from the surface, 


(4): -(6) 2 


y= Ture u\?’ 
a), ~ (a), 


where (u/U)p is the reading of the Pitot tube. 

A suitable procedure for determining the y-distribution by this method 
is to traverse a Pitot tube along the transition region at a fixed height above 
the surface. ‘To get the best results, the height should be such that the 
difference in the velocity, at that height, between the laminar and turbulent 
profiles is as large as possible throughout the transition region. From an 
examination of figures 9 and 10, it seems that this height should be of the 
order of 6. A surface tube is not likely to give satisfactory results, as its 
readings cannot always be quantitatively interpreted. 


y can be obtained as 
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Figure 16 compares intermittency measurements obtained for the 
same flow by the hot wire and Pitot methods. For the latter, the velocities 
(u/U), were read off, for a height from the surface of 0-024 in., from 
measured Pitot velocity profiles. Figure 17 shows the results obtained 
from a Pitot traverse at a height of 0-035 in., together with the calculated 
laminar and turbulent curves. The origin for the turbulent boundary 
layer was chosen as the point where the measured curve departs from the 
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Figure 16. Estimation of y using a Pitot tube. 
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Figure 17. Pitot traverse along plate to estimate } 


laminar. Incidentally, a simple method of obtaining the ‘turbulent’ curve 
in figure 17 is to plot the measured fully-turbulent points on logarithmic 
paper against (v—.«,) and draw a straight line through them extending 
upstream. This line then gives the (u/U) , values for the y determination 
in the transition region. ‘This procedure gives good results and would be 
completely justified if the turbulent profiles followed a power law with a 
constant index throughout. The values of the intermittency obtained 
from this experiment are also plotted in figure 16. 
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CONCLUDING REMARKS 
In spite of the origin of turbulence remaining obscure, Emmons’ picture 
of transition in a laminar boundary layer supplemented by experimental 
information allows a detailed physical model of the process to be set up. 
From this a relatively adequate calculation of the macroscopic effects can 
be carried out. 


One of us (R. N.) wishes to acknowledge his indebtedness to the National 
Institute of Sciences of India for a Fellowship held during the course of 
this research. 
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REVIEWS 


Jets, Wakes and Cavities, by Garretr BirkuorF and E. H. ZARANTONELLO. 
New York: Academic Press, 1957. 353 pp. $10.00 or 80s. 


There can be little doubt that Helmholtz, with his paper ‘‘ Uber 
discontinuerliche Fliissigkeitsbewegungen ”’ (1868), laid one of the corner- 
stones of modern fluid mechanics. He himself, with the editors of The 
Philosophical Magazine, who published a translation of the paper in the 
same year, was confident of its importance—just as, thirty-six years later, 
Prandtl was confident in his paper on the boundary layer. Yet, Helmholtz 
seems not to have known which cornerstone he was laying, and probably 
he imagined that he was laying Prandtl’s. 

Certainly the discontinuous solutions of Laplace’s equation constructed 
by Helmholtz, and by his successors such as Kirchhoff (who introduced the 
name ‘ free streamline ’) and Rayleigh, have found effective application almost 
exclusively to flows of liquid bordered by regions of gas or vapour—notably, 
liquid jets in a gaseous atmosphere, and liquid flows about obstacles with 
cavity formation—whereas the flows of a single homogeneous fluid, to 
which the authors originally intended the solutions to apply, have proved too 
complicated for detailed study by these means. Meanwhile, the theory 
and practice of constructing such discontinuous solutions have been 
developed to a high degree of virtuosity by an army of workers, and many 
interesting phenomena have been illuminated as a result. Nevertheless, 
too few of those writing on the subject have grasped its true relation to 
physical phenomena, and many have repeated almost word for word the 
errors of the founders of the theory, such as Helmholtz’s ascription of 
flow separation from a solid surface to causes which really are those under- 
lying cavitation. Especial harm has been done by the writers of hydro- 
dynamical textbooks, who have actively propagated false views of the 
significance of free-streamline theory. 

The vast structure built upon Helmholtz’s original paper—comprising 
the mathematical theory of the existence and uniqueness of flows with free 
streamlines ; the techniques of obtaining analytical expressions for such 
flows in particular cases and then translating them into numerical results ; 
the use of the solutions to throw light on jet collisions, on planing, on the 
effectiveness of shaped charges, on vapour-filled and air-filled cavities, 
on contraction coefficients and flows over weirs and waterfalls ; the 
experimental and theoretical study of the physical conditions in flows 
involving jets and cavities and the relation between these conditions and 
aspects of the mathematical solutions ; the further study of possible 
relationships between separated flows of a homogeneous fluid and particular 
types of free-streamline solution—all this is a structure spanning the widely 
different disciplines of pure and applied mathematics, and involving 
substantial parts of physics and engineering. ‘The need for a unified 
account of all this varied material, which would help those working on 
isolated parts of it to see the material as a whole, and reduce the danger of 
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their misinterpreting its significance, has long existed, as has the difficulty of 
finding anyone competent to undertake such a synthesis. 

Now Professors Garrett Birkhoftf and E. H. Zarantonello, both 
distinguished workers in the field, have collaborated to produce just such a 
unified account, and it is a real pleasure to report that they have succeeded, 
achieving both the necessary breadth and the necessary depth. Professor 
Birkhoff is to be congratulated on having made the difficult transition from 
pure to applied mathematics, at least for the purpose of writing parts of 
this book; and those applied mathematicians who found their enjoyment of 
his Hydrodynamics spoilt by its oversimplified view of the relations 
between natural science and mathematics can be assured that this experience 
will not be repeated. 

The authors are particularly strong adherents to the view that free 
streamlines are valuable only as representations of liquid-gas boundaries. 
With the principal object of emphasizing this, they include three chapters 
on wakes and jets in the flow of homogeneous fluid (for example, air jets in 
air), chapters which summarize rather briefly the knowledge in these fields 
but at least make clear to any doubters that an approach to them by means 
of free-streamline solutions must be quite inadequate. With the exception 
of these chapters (12-14) the book is a study of cavities and liquid jets, and of 
the theory of the mathematical determination of free-streamline solutions. 

The authors are careful not to give an oversimplified view of the cavitation 
phenomenon, and they devote their final chapter to its many complexities 
as well as to descriptions of the various instabilities which lead to the 
break-up of jets. 

Naturally, one of the most useful features of the book is the ‘ zoo’ of 
basic free-streamline flows, that is, two-dimensional irrotational flows with 
either the magnitude or the direction of the velocity uniform along each 
portion of the boundary. ‘The zoo is well arranged, in three chapters 
(2, 3 and 5) dealing with types of flow in increasing order of complexity. 
In these chapters there are some sixty sketches of particular configurations, 
showing simply the solid boundaries and the free streamlines in each case, 
which, with a sentence or two in the text, identify the flow clearly. In 
addition, the analytical expressions for the solution are given, and remarks 
about its utility, if any. In chapter 9 the numerical computation of these 
flows is systematically discussed. Streamline shapes and other details of 
the internal flow are obtained in many cases and exhibited in a series of 
plates at the back of the book. ‘This is a big advance, since, in previous 
work, details only of the boundary streamlines had been computed for flows 
of this kind. In addition, many tables and graphs are given of the relations 
between the defining parameters for particular types of flow. 

The general theory of the existence, uniqueness and nature of free- 
streamline flows, including those with curved solid boundaries, is treated in 
chapters 4,6 and 7. A knowledge of at least the results of this material is of 
the greatest help to those using free-streamline solutions. Readers 
wishing to understand the proofs will have to make a very real effort, 
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but fortunately they are set out in a civilized way which does not prevent a 
reader skimming the cream of the results while only sipping at the difficult 
ideas underlying their demonstration. 

The first half of chapter 8 is devoted to the Chaplygin compressible- 
flow free-streamline theory of gaseous jets, a subject included somewhat 
inconsistently in view of the authors’ general contention that only liquid 
jets are well represented by free-streamline theory. ‘This chapter ends with 


some remarks about the problem of including gravitational effects. 
However, it gives no account of perturbation methods in this connection, 
and computations by the relaxation method are not mentioned, being only 
briefly treated in a section at the end of chapter 9, whereas they have 
proved especially fruitful for calculating free streamline flows under 
gravity. Chapter 10 is a good account of what is known about axisymmetric 
free-streamline flows, and chapter 11 takes us rapidly, but with references, 
through an extensive miscellany of unsteady flows with free boundaries. 

When one looks back over this book, which is such a welcome addition 
to the shelf of really good hydrodynamical treatises, it is difficult to find any 
general questions on which one would like to challenge the authors. How- 
ever, I willtry. Perhaps the one possible fault is an inevitable consequence 
of one of the book’s chief virtues, namely its insistence (so necessary in 
view of the errors of the past) that flows of homogeneous fluid with 
separation are too complicated to be represented by free-streamline solutions. 
At the end of the book, on the other hand, one would like to ask: ‘‘ Given 
that this is indisputably true in the main, are there not still occasional 
points where free-streamline theory can throw some light on flows of 
homogeneous fluid ?”” In other words, was Helmholtz completely wrong ¢ 

For example, the authors’ figure 3 in chapter 2 indicates that the pressure 
distribution on a flat plate, at an angle of incidence of 15° or more, is close 
to that given by one of the cavitation models (such as that with a re-entrant 
jet, but admittedly all the models give much the same distribution over the 
plate), with a cavity pressure equal to that found experimentally on the 
leeward side of the plate. Also, figures 2 in both chapter 1 and chapter 14 
indicate that a re-entrant jet model is a reasonable qualitative representation 
of the flow near the plate, although the real flow becomes fuzzier and fuzzier 
(relative to the model) as one goes away from the plate. 

Again, there are experimental indications that vortex sheets in the 
neighbourhood of solid boundaries are fairly stable, so that free-streamline 
models may be more useful for representing ‘ separation bubbles ’ (that is, 
finite dead-air regions produced by boundary-layer separation followed by 
reattachment) than, say, wakes. ‘The choice of model for a separation 
bubble would involve boundary-layer considerations to determine the 
position of separation. 

But it probably needs to be pointed out that the authors of this book are 
at their weakést in boundary-layer considerations. Fortunately, the main 
sections of the book do not lose much from this, since on most solid 
boundaries in the free-streamline solutions which are discussed the flow is 
accelerating, so that boundary layers should be thin and should not separate. 
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However, there are several flows illustrated involving flow round 
a corner in a wall (see chapter 2, figures 4a, 4b, 6d; chapter 3, 
figures 10b, 12a ; chapter 5, figures 7a, 7b, 8a, 8b, 8c, 8d, 8e, 10a, 10b, 
10c). In all of these, boundary-layer separation would cause large 
inaccuracies in the solution—a circumstance to which the reader’s 
attention should have been drawn more explicitly. 

Actually, there is no evidence in the book itself that the authors are 
aware of this fact in the case of corners concave to the flow. If they had 
been, they would presumably not have delayed the only discussion of 
boundary layers and separation until chapter 14, where a single page is 
devoted to the subject. This page, indeed, is a peculiar one, apparently 
designed to play down the role of the boundary layer, and it fails to state 
clearly that separation will not occur at all in flows with sufficiently 
small adverse pressure gradients along solid boundaries. 

It is surprising that authors who have written at such length on irrota- 
tional flows should wish to play down the boundary layer, by whose presence 
alone irrotational flows are made possible in the face of the hydrodynamic 
boundary condition of zero relative velocity at a solid boundary. ‘The 
Navier-Stokes equations for incompressible flow are exactly satisfied by irro- 
tational flow fields, but these are compatible with the boundary condition 
only if a boundary layer intervenes between them and any solid boundary. 
A flow irrotational outside this thin layer will be possible only if the equa- 
tions governing the dynamics of the layer, in the environment of the given 
irrotational flow, permit a steady solution without such a singularity as is 
implied by the separation condition. Divergences between theory and 
experiment for boundary layers which do separate (that is, for which the 
equations reach this singularity) constitute no criticism of a theory whose 
most important conclusion is the condition under which separation will be 
absent. 

However, this condition is not satisfied by any flow along a wall up to a 
stagnation point (like a concave corner). It is for this reason that flows 
involving separation bubbles (see for example chapter 5, figures 2b, 2c, 2d in 
this book) have been proposed to describe such cases. ‘The present 
reviewer has argued for such a model in the case of flow up a step (appar- 
ently without convincing Professor Birkhoff, for one cannot, alas, suppose 
the reference unknown to such an accomplished bibliographer!). In 
addition, there is reason to hope that free-streamline theories of separa- 
tion bubbles on thin aerofoils at high angle of incidence may gradually 
achieve greater accuracy. 

Still, it must be confessed that the authors’ case remains a tenable one, 
and that Helmholtz’s vindication is yet to come. In the meantime, the 
appearance of this powerful critical survey of the present position in the 
whole subject may be regarded as bearing witness to the vitality of his 
original conception. Our thanks are due to the authors for making such a 
book as pleasant to read as it is comprehensive and instructive. 


M. J. LIGHTHILL 
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